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ABSTRACT 


Experiments  were  performed  to  investigate  condensation  of  FC-72  along  parallel, 
square  micro-channels  with  a  hydraulic  diameter  of  1  mm  and  a  length  of  29.9  cm,  which 
were  formed  in  the  top  surface  of  a  solid  copper  plate.  The  condensation  was  achieved 
by  rejecting  heat  to  a  counter  flow  of  water  through  channels  brazed  to  the  underside  of 
the  copper  plate.  The  FC-72  entered  the  micro-channels  slightly  superheated,  and 
operating  conditions  included  FC-72  mass  velocities  of  68  -  367  kg/m2s,  FC-72  saturation 
temperatures  of  57.2  -  62.3 °C,  and  water  mass  flow  rates  of  3  -  6  g/s.  Using  high-speed 
video  imaging  and  photomicrographic  techniques,  five  distinct  flow  regimes  were 
identified:  smooth-annular,  wavy-annular,  transition,  slug,  and  bubbly,  with  the  smooth- 
annular  and  wavy-annular  regimes  being  most  prevalent. 

A  detailed  pressure  model  is  presented  which  includes  all  components  of  pressure 
drop  across  the  micro-channel.  Different  sub-models  for  the  frictional  and  accelerational 
pressure  gradients  are  examined  using  the  homogenous  equilibrium  model  (with  different 
two-phase  friction  factor  relations)  as  well  as  previous  macro-channel  and  mini/micro- 
channel  separated  flow  correlations.  Unexpectedly,  the  homogenous  flow  model 
provided  far  more  accurate  predictions  of  pressure  drop  than  the  separated  flow  models. 
Among  the  separated  flow  models,  better  predictions  were  achieved  with  those  for 
adiabatic  and  mini/micro-channels  than  those  for  flow  boiling  and  macro-channels. 

The  FC-72  condensation  heat  transfer  coefficient  was  highest  near  the  channel 
inlet,  where  the  annual  liquid  film  is  thinnest.  The  heat  transfer  coefficient  decreased 
along  the  micro-channel  because  of  the  film  thickening  and  eventual  collapse  of  the 
annular  regime.  Notable  heat  transfer  enhancement  was  observed  for  annular  flow 
regions  of  the  micro-channel  associated  with  interfacial  waves.  Comparing  the  present 

data  to  predictions  of  previous  annular  condensation  heat  transfer  correlations  shows 
correlations  intended  for  macro-channels  generally  provide  better  predictions  than 
correlations  intended  specifically  for  mini/micro-channels.  A  new  condensation  heat 
transfer  coefficient  correlation  is  proposed  for  annular  condensation  heat  transfer  in 
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mini/micro-channels.  The  new  correlation  shows  excellent  predictive  capability  based  on 
both  the  present  FC-72  data  and  a  large  database  for  mini/micro-channel  flows  amassed 
from  eight  previous  sources. 

A  theoretical  control-volume-based  model  is  proposed  based  on  the  assumptions 
of  smooth  interface  between  the  annular  liquid  film  and  vapor  core,  and  uniform  film 
thickness  around  the  channel’s  circumference.  Mass  and  momentum  conservation  are 
applied  to  control  volumes  encompassing  the  liquid  film  and  the  vapor  core  separately. 
The  model  accounts  for  interfacial  suppression  of  turbulent  eddies  due  to  surface  tension 
with  the  aid  of  a  new  eddy  diffusivity  model  specifically  tailored  to  shear-driven 
turbulent  films.  The  model  predictions  are  compared  with  present  experimental  pressure 
drop  and  heat  transfer  data  for  annular  condensation  of  FC-72  mass  velocities  of  248  - 
367  kg/m2s,  saturation  temperatures  of  57.8  -  62.3°C,  qualities  of  0.23  -  1.0,  and  water 
mass  flow  rates  of  3  -  6  g/s.  The  data  are  also  compared  to  predictions  of  previous 
separated  flow  mini/micro-channel  and  macro-channel  correlations.  While  some  of  the 
previous  correlations  do  provide  good  predictions  of  the  average  heat  transfer  coefficient, 
they  fail  to  capture  axial  variation  of  the  local  heat  transfer  coefficient  along  the  channel. 
The  new  model  accurately  captures  the  pressure  drop  and  heat  transfer  coefficient  data  in 
both  magnitude  and  trend,  evidenced  by  mean  absolute  error  values  of  3.6%  and  9.3%, 
respectively. 
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CHAPTER  1.  INTRODUCTION 


1.1.  Motivation  for  Study 

Developments  in  many  cutting-edge  technologies  are  becoming  increasingly 
dependent  upon  the  ability  to  dissipate  large  amounts  of  heat  from  small  surface  areas. 
Examples  include  high  performance  computers,  hybrid  vehicle  power  electronics,  lasers, 
radars  and  avionics.  Although  single-phase  cooling  systems  have  been  successfully 
implemented  in  the  past,  the  fast  increase  in  power  density  in  these  technologies  is 
pushing  single-phase  systems  into  unchartered  territory  in  terms  coolant  flow  rate  and 
pressure  drop,  let  alone  the  need  to  greatly  increase  heat  transfer  area.  Even  with  these 
provisions,  cooling  goals  in  many  applications  are  becoming  virtually  impossible  to 
achieve  with  single-phase  systems.  Because  of  these  inherent  limitations,  focus  has 
shifted  in  recent  years  in  favor  of  two-phase  cooling  schemes. 

Modern  technologies  are  demanding  more  effective  schemes  to  tackle  heat 
removal  from  very  high  power  density  devices.  Most  of  the  phase-change  electronic 
cooling  research  published  during  the  past  three  decades  has  been  focused  on  removing 
heat  from  the  electronic  device  using  micro-channel  heat  sinks,  jet-impingement  or 
sprays  (Mudawar  and  Wadsworth,  1991;  Jiang  et  al.,  2001;  Kawahara  et  al.  2002;  Lin 
and  Ponnappan,  2003;  Lee  and  Mudawar,  2005).  Far  less  emphasis  has  been  placed  on 
high-flux  heat  rejection  from  the  two-phase  cooling  system.  A  key  reason  behind  this 
trend  is  a  common  perception  that  a  commercial  condenser  can  always  be  found  to  reject 
the  heat  from  virtually  any  phase-change  cooling  system.  However,  recent  studies 
concerning  very  high  power  density  defense  electronics  have  shown  that  substantial 
enhancement  in  heat  dissipation  from  the  device  can  be  achieved  by  using  a  two-loop 
cooling  system  (Lee  and  Mudawar,  2008a,  2008b).  Here,  a  primary  cooling  loop  is  used 
to  extract  the  heat  from  the  device  by  highly  subcooled  flow  boiling,  and  the  heat  is 
rejected  via  a  high  performance  heat  exchanger  to  a  separate  low-temperature 
refrigeration  loop.  Compact  and  lightweight  system  design  requires  that  the  intermediate 
condenser  achieve  heat  fluxes  comparable  to  those  of  the  flow  boiling  module. 
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1.2.  Miniature  Condenser  for  High-Flux  Heat  Rejection 

The  superior  cooling  performance  of  two-phase  cooling  systems  is  realized  by  the 
enormous  boiling  and  condensation  heat  transfer  coefficients  associated  with  heat 
acquisition  from  the  heat-dissipating  device  and  heat  rejection,  respectively.  The  quest 
for  compact  and  lightweight  packaging  favors  the  use  of  miniature  boiler  designs,  such  as 
micro-channel  heat  sinks  (Lee  and  Mudawar,  2008a,  2008b)  or  jet-impingement 
(Mudawar  and  Wadsworth,  1991).  While  many  present  two-phase  cooling  systems 
employ  fairly  standard  air-cooled  condensers  to  reject  the  heat  to  the  ambient,  there  is 
now  a  growing  need  for  miniature  condensers  that  can  reject  the  heat  by  condensing  a 
primary  coolant  in  a  compact  primary  cooling  loop.  The  heat  is  transferred  to  a 
secondary  liquid  coolant  and  transported  to  a  remote  heat  exchanger  where  is  ultimately 
rejected  to  ambient  air  (or  seawater  for  marine  applications).  As  described  by  Lee  and 
Mudawar  (2008a,  2008b),  the  secondary  cooling  loop  could  consist  of  a  low-temperature 
refrigeration  cycle.  Using  separate  cooling  loops  decouples  the  performance  of  the  boiler 
from  that  of  the  refrigeration  cycle,  enabling  the  dissipation  of  much  higher  device  heat 
fluxes  than  with  an  evaporator  incorporated  directly  in  the  refrigeration  loop.  Integrating 
miniature  boilers  and  condensers  in  the  compact  primary  cooling  loop  will  necessitate 
greatly  enhancing  the  condenser’s  performance  to  a  level  commensurate  with  that  of  the 
boiler.  This  goal  can  be  realized  with  a  condenser  containing  parallel  micro-channels 
similar  to  those  utilized  in  micro-channel  flow  boiling  heat  sinks.  However,  special  care 
must  be  exercised  to  achieve  the  thermal  goals  without  compromising  the  system’s 
efficiency  by  the  increased  pressure  drop. 

Because  condensation  heat  transfer  coefficients  are  typically  much  smaller  than 
those  realized  with  subcooled  flow  boiling,  commercial  condensers  are  often  far  too  large 
to  meet  the  size  and  packaging  constraints  of  defense  electronics.  The  present  study 
concerns  the  design  of  a  relatively  new  class  of  condensers  that  employ  a  series  of  micro¬ 
channels  to  meet  the  stringent  size  and  weight  requirements  of  defense  electronics. 
Another  goal  is  to  maintain  mostly  annular  flow  along  the  micro-channels  to  capitalize 
upon  the  large  condensation  heat  transfer  coefficients  associated  with  thin  films.  With 
superheated  or  saturated  inlet  conditions,  a  very  thin  film  is  initiated  in  the  upstream 
region  of  the  channel,  which  is  driven  along  the  channel  by  the  shear  stresses  exerted  by 
the  core  vapor  flow.  Micro-channels  greatly  increase  vapor  velocity  and  therefore  the 


Boiling  and  Two-Phase  Flow  Laboratory 


3 


shear  stress  exerted  upon  the  film  interface.  This  greatly  decreases  the  film  thickness, 
resulting  in  very  high  condensation  heat  transfer  coefficients.  However,  the  film 
thickness  increases  along  the  micro-channel  as  an  increasing  fraction  of  the  vapor  flow 
condenses  into  liquid.  Eventually,  the  annular  regime  collapses  and  is  replaced  by  a 
succession  of  thermally  less  efficient  slug,  bubbly,  and  liquid  flow  regimes.  While 
micro-channels  do  enhance  heat  transfer  performance,  they  also  increase  pressure  drop. 
Therefore,  the  design  of  micro-channel  condensers  requires  predictive  tools  for  both 
pressure  drop  and  condensation  heat  transfer  coefficient. 


1.3.  Transport  Phenomena  in  Micro-Channel  Condenser 

Two-phase  pressure  drop  is  comprised  of  frictional  as  well  as  accelerational  or 
decelerational  components;  acceleration  is  associated  with  boiling  flows  and  deceleration 
condensing  flows.  In  high-flux  micro-channel  boiling  flows,  acceleration  contributes  a 
large  fraction  of  the  total  pressure  drop  and,  therefore,  pumping  power.  However,  in 
high-flux  micro-channel  condensing  flows,  deceleration  provides  the  benefit  of 
decreasing  the  pressure  drop. 

Another  factor  that  influences  pressure  drop  in  condensing  flows  is  two-phase 
flow  regime.  Condensing  flow  regimes  in  a  tube  can  generally  be  classified  into  pure 
vapor,  annular,  slug,  bubbly  and  pure  liquid.  The  annular  flow  regime  is  the  most  crucial 
for  condenser  design  since,  for  most  practical  designs,  it  prevails  over  the  largest  fraction 
of  the  tube  length  and  provides  very  high  heat  transfer  rates.  The  annular  flow  regime 
commences  upstream  in  the  tube  once  vapor  begins  to  condense  into  a  liquid  film  along 
the  tube  wall.  This  film  is  typically  very  thin  and  associated  with  a  small  thermal 
resistance  between  the  vapor  and  the  wall,  resulting  in  very  high  heat  transfer  coefficients 
compared  to  the  other  flow  regimes.  As  more  of  the  vapor  is  condensed  into  liquid,  the 
annular  structure  collapses  and  the  slug  and  bubbly  regimes  are  formed  in  succession. 

Key  to  achieving  high  condensation  heat  transfer  coefficients  is  to  maintain 
annular  flow  over  a  large  fraction  of  the  micro-channel's  length.  This  flow  regime 
consists  of  a  liquid  film  that  sheathes  the  channel  walls,  driven  by  the  shear  force  exerted 
by  the  vapor  core.  The  film  is  extremely  thin  in  the  upstream  region  of  the  micro-channel 
where  the  annular  regime  is  initiated,  which  yields  very  high  condensation  heat  transfer 
coefficients.  The  heat  transfer  coefficient  decreases  along  the  stream-wise  direction  as 
the  annular  film  gradually  thickens  and,  especially,  as  the  annular  regime  is  replaced,  in 
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succession,  by  transition,  slug,  bubbly  and  single-phase  liquid  flow  regimes. 
Recognizing  the  importance  of  annular  flow,  previous  condensation  studies  have  focused 
mostly  on  this  particular  regime  (e.g.,  Akers  et  al.,  1958;  Cavallini  and  Zecchin,  1974; 
Shah,  1979;  Dobson  and  Chato,  1998;  Moser  et  al.,  1998;  Want  et  al.,  2002;  Koyama  et 
al.,  2003;  Huang  et  al.,  2010). 

The  pursuit  of  compact  condenser  design  is  evident  from  several  recent  articles 
concerning  multi-port  micro-channel  condensers  (Yang  and  Webb,  1996;  Yan  and  Lin, 
1999;  Wang  et  al.,  2002;  Koyama  et  al.,  2003;  Park  and  Hrnjak,  2009).  Koyama  et  al. 
(2003)  examined  the  performance  of  a  condenser  having  rectangular  micro-channels  with 
a  hydraulic  diameter  of  0.81  mm  using  R134a  as  working  fluid.  The  highest  heat  transfer 
coefficient,  on  the  order  of  104  W/m2K,  was  achieved  in  the  high  quality  region  with  a 
relatively  high  mass  velocity  of  G  =  652  kg/m2s.  Recently,  Park  and  Hrnjak  (2009) 
examined  the  condensation  of  CO2  in  0.89-mm  circular  parallel  micro-channels  and 
achieved  maximum  heat  transfer  coefficients  on  the  same  order  as  Koyama  et  al.  with  G 
=  600  kg/m2s. 

While  past  studies  have  provided  valuable  databases  for  condensation  in  micro¬ 
channels,  there  is  a  general  lack  of  fundamental  understanding  of  interfacial  behavior  and 
shortage  of  predictive  tools  for  both  pressure  drop  and  heat  transfer  coefficient.  On  the 
fundamental  side,  two  key  phenomena  that  are  of  great  interest  to  modeling  condensation 
in  micro-channels  are  interfacial  waves  and  turbulence  in  the  annular  liquid  film.  Past 
studies  involving  adiabatic,  heated  and  evaporating  liquid  films  have  shown  interfacial 
waves  can  greatly  influence  mass,  momentum  and  heat  transfer  in  the  film  (Lyu  and 
Mudawar,  1991a,  1991b;  Mudawar  and  Houpt,  1993a,  1993b;  Shmerler  and  Mudawar, 
1998a,  1998b).  Surface  tension  at  the  vapor-liquid  interface  has  been  shown  by 
Mudawar  and  El-Masri  (1986)  to  cause  significant  dampening  of  turbulent  fluctuations  in 
films,  which  can  have  a  profound  influence  on  heat  transfer  across  the  film.  Aside  from 
these  fundamental  issues,  there  is  also  a  need  for  reliable  correlations  for  condensing 
micro-channel  flows.  These  issues  and  needs  are  the  primary  motivation  for  the  present 
study. 


1.4.  Two-Phase  Flow  Model  for  Annular  Flow  Condensation 

The  vast  majority  of  published  studies  addressing  the  prediction  of  pressure  drop 
and  condensation  heat  transfer  coefficient  for  annular  flows  are  based  on  semi-empirical 
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separated  flow  correlations  (mostly  Lockhart-Martinelli-type,  1949)  and  empirical  heat 
transfer  correlations  ( e.g Shah,  1979),  respectively.  Like  all  correlations,  these 
predictive  tools  are  valid  only  for  the  range  of  operating  conditions  of  databases  these 
correlations  are  based  upon.  This  limitation  has  created  a  need  for  development  of 
theoretical  models  that  possess  a  broader  application  range.  Theoretical  control-volume- 
based  models  have  been  quite  effective  in  predicting  both  pressure  drop  (Qu  and 
Mudawar,  2003a)  and  heat  transfer  (Qu  and  Mudawar,  2003c)  in  saturated  flow  boiling  in 
micro-channel  heat  sinks,  subcooled  flow  boiling  pressure  drop  in  micro-channel  heat 
sinks  (Lee  and  Mudawar,  2009),  and  flow  boiling  critical  heat  flux  (Sturgis  and  Mudawar, 
1999).  Nonetheless,  several  fundamental  challenges  remain  when  attempting  to  develop 
an  accurate  model  for  annular  two-phase  flow.  These  include  interfacial  instabilities, 
interfacial  mass,  momentum  and  heat  transfer,  and  turbulence  within  the  annular  film. 

One  source  of  difficulty  in  modeling  annular  two-phase  flow  is  interfacial  waves 
(Hubbard  et  at.,  1976;  Ueda  et  a/.,  1977;  Shmerler  and  Mudawar,  1988a,  1988b;  Koskie 
et  ai,  1989;  Lyu  and  Mudawar,  1991a;  Mudawar  and  Houpt,  1993b).  These  are  highly 
complex  phenomena  and  no  effective  method  has  been  developed  for  their 
characterization.  One  type  of  waves  is  ripples,  which  are  characterized  by  both  small 
amplitude  and  small  wavelength  and  do  not  contribute  substantially  to  liquid  film  mass 
transport.  Large  waves  can  also  develop,  where  the  perturbed  liquid  film  is  replaced  by 
lumps  of  liquid  that  are  separated  by  a  thin  layer  of  liquid.  Due  to  their  appreciable 
protrusion  into  the  vapor  flow,  the  motion  of  large  waves  can  be  dominated  more  by 
vapor  drag  forces  than  by  interfacial  vapor  shear.  Modeling  large  waves  is  complicated 
by  the  fact  that  they  are  highly  nonlinear  forms  of  instability  that  greatly  complicate  the 
effectiveness  of  averaging  methods  to  predict  the  liquid  film’s  momentum  or  heat  transfer 
transport  behavior.  Despite  extensive  efforts  to  assess  the  influence  of  interfacial  waves 
on  mass,  momentum  and  heat  transfer,  no  systematic  tools  have  been  developed  to 
effectively  model  this  influence. 

Two  other  phenomena  that  complicate  the  development  of  accurate  annular  flow 
models  are  droplet  entrainment  and  droplet  deposition.  These  two  phenomena  are  closely 
related  to  interfacial  waviness.  Entrainment  refers  to  tiny  liquid  droplets  that  break  off 
the  crests  of  waves  and  are  driven  along  with  the  vapor  core.  Deposition  refers  to 
droplets  from  the  vapor  core  falling  back  upon  the  film  interface.  In  their  study  of 
annular  film  evaporation,  Qu  and  Mudawar  (2003a)  proved  that  droplet  entrainment  and 
deposition  play  a  very  important  role  in  micro-channels  compared  to  macro-channels. 
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Another  challenge  in  modeling  annular  flows  is  accurate  prediction  of  turbulent  or 
eddy  diffusivity  within  the  annular  liquid  film.  Classical  turbulence  models  are  known  to 
break  apart  when  dealing  with  fluid  regions  near  a  vapor-liquid  interface  (Hubbard  et  al., 
1976;  Ueda  et  al.,  1977;  Mudawar  and  El-Masri,  1986).  As  discussed  by  Mudawar  and 
El-Masri,  surface  tension  forces  along  the  interface  can  significantly  dampen  turbulent 
eddies,  resulting  in  appreciable  resistance  to  heat  transfer  near  the  interface.  Ignoring  this 
effect  could  result  in  appreciable  error  in  predicting  the  condensation  heat  transfer 
coefficient. 


1.5.  Research  Objectives 

This  study  examines  the  two-phase  flow  regimes,  pressure  drop  and  heat  transfer 
associated  with  condensation  of  FC-72  along  parallel  micro-channels  with  a  hydraulic 
diameter  of  1  mm  (Kim  et  al.,  in  review;  Kim  and  Mudawar,  in  review-a).  The 
construction  of  the  condensation  module,  two-phase  condensation  loop  and  experimental 
methods  used  will  be  described.  Using  high-speed  video  imaging  and  photomicrographic 
techniques,  five  distinct  flow  regimes  were  identified:  smooth-annular,  wavy-annular, 
transition,  slug,  and  bubbly,  with  the  smooth-annular  and  wavy-annular  regimes  being 
most  prevalent,  which  are  compared  to  those  from  previous  studies.  The  FC-72  pressure 
drop  data  are  compared  with  the  predictions  of  prior  two-phase  homogenous  (Bottomley, 
1936;  Lewis  and  Robertson,  1940;  Benjamin  and  Miller,  1942;  Markson  et  al.,  1942; 
McAdams  et  al.,  1942;  Allen,  1951;  Akers  et  al.,  1958;  Cicchitti  et  al,  1960;  Owens, 
1961;  Dukler  et  al.,  1964;  Beattie  and  Whalley,  1982;  Lin  et  al.,  1991)  and  separated 
flow  models  (Lockhart  and  Martinelli,  1949;  Friedel,  1979;  Chisholm,  1983;  Muller- 
Steinhagen  and  Heck,  1986;  Jung  and  Radermacher,  1989;  Mishima  and  Hibiki,  1996; 
Wang  et  al.,  1997;  Tran  et  al.,  2000;  Chen  et  al.,  2001;  Lee  and  Lee,  2001;  Zhang  and 
Webb,  2001;  Cavallini  et  al.,  2002;  Yu  et  al.,  2002;  Lee  and  Mudawar,  2005;  Nino  et  al., 
2005;  Hwang  and  Kim,  2006)  to  identify  most  suitable  models  for  micro-channel 
condensing  flows.  The  present  FC-72  heat  transfer  data  are  first  compared  with 
predictions  of  several  previous  correlations  for  both  macro-channels  (Akers  et  al.,  1958; 
Cavallini  and  Zecchin,  1974;  Shah,  1979;  Dobson  and  Chato,  1998;  Moser  and  Webb, 
1998)  and  mini/micro-channels  (Wang  et  al.,  2002;  Koyama  et  al.,  2003;  Huang  et  al., 
2010).  A  new  correlation  for  annular  condensation  heat  transfer  is  proposed,  which  is 
validated  against  the  present  data  in  addition  to  a  database  for  condensation  in 
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mini/micro-channels  from  eight  previous  sources  (Hirofumi  and  Webb,  1995;  Zhang, 
1998;  Yan  and  Lin,  1999;  Wang  et  al,  2002;  Cavallini  et  al.,  2005;  Matkovic  et  al., 
2009;  Park  and  Hrnjak,  2009;  Huang  et  al .,  2010). 

Another  goal  of  the  present  study  is  to  develop  a  theoretical  control-volume-based 
model  for  annular  condensing  micro-channel  flows  (Kim  and  Mudawar,  in  review-b). 
Addressed  in  the  development  of  this  new  model  are  the  effects  of  droplet  entrainment 
and  deposition,  interfacial  instabilities,  and  dampening  of  liquid  film  turbulence  near  the 
liquid-vapor  interface.  The  model’s  predictive  accuracy  is  assessed  against  the  FC-72 
pressure  drop  and  heat  transfer  data  measured  previously  by  the  authors  (Kim  et  al.,  in 
review;  Kim  and  Mudawar,  in  review-a),  and  compared  to  those  of  prior  annular 
condensation  correlations. 

Another  goal  is  to  extend  both  the  theoretical  models  and  correlation  approaches 
to  flow  boiling  in  micro-channels. 
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CHAPTER  2.  EXPERIMENTAL  METHODS 

2.1.  Condensation  Flow  Loop 

Figure  2.1  shows  a  schematic  of  the  two-phase  experimental  apparatus 
constructed  for  this  study.  It  uses  two  sub-loops,  one  for  FC-72,  the  primary  coolant,  and 
the  second  for  the  cooling  water.  A  3M-company  Fluorinet  electronic  liquid,  FC-72  has 
excellent  dielectric  properties,  a  relatively  moderate  boiling  point  of  56°C  at  1  bar,  and  a 
surface  tension  much  smaller  than  that  of  water,  0.0084  compared  to  0.059  N/m  at  1  bar. 
Table  2.1  provides  a  summary  of  the  thermophysical  properties  of  FC-72  at  Tsat  =  60°C, 
which  the  saturation  temperature  corresponding  an  operating  pressure  close  to  that  used 
in  the  present  experiments. 

The  FC-72  is  circulated  through  the  primary  loop  with  the  aid  of  a  gear  pump. 
The  fluid  passes  through  a  set  of  flow  meters  followed  by  an  in-line  electric  heater,  which 
is  powered  by  a  variable  voltage  transformer,  before  entering  the  micro-channel 
condensation  module.  Both  temperature  and  pressure  are  measured  in  the  inlet  plenum  of 
condensation  module  to  determine  fluid  quality,  which  can  be  regulated  to  the  desired 
value  with  the  aid  of  a  variable  voltage  transformer  that  powers  the  in-line  heater 
upstream.  The  FC-72  condenses  along  the  condensation  module  by  rejecting  heat  to  the 
waterside.  The  FC-72  temperature  and  pressure  are  measured  once  more  at  the 
condensation  module’s  outlet  plenum.  Exiting  the  condenser  module,  the  FC-72  is 
passed  through  a  plate-type  heat  exchanger  to  condense  any  remaining  vapor  and  achieve 
the  desired  temperature  as  the  FC-72  returns  to  the  primary  loop’s  reservoir. 


Table  2.1  Thermophysical  properties  of  FC-72  at  Tsa,  =  60°C. 
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Pg 

Mf 
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cpJ 
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kJ/kg 

kg/m3 

kg/m3 
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W/rn-K 

kJ/kg-K 

mN/m 

93.7 

1583.4 

14.90 

4. 18x1  O'4 
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Figure  2.1  Schematic  diagram  of  test  loop. 
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In  the  condensation  module,  water  travels  in  a  counter-flow  direction  to  that  of  the 
FC-72.  The  water  enters  the  condensation  module  at  near  room  temperature  and  warms 
up  by  capturing  heat  from  the  condensing  FC-72  flow.  The  warm  water  flows  through  a 
Lytron  modular  cooling  system,  where  it  cools  back  to  near  room  temperature.  The  water 
flow  exiting  the  Lytron  cooling  system  is  divided  into  two  parts,  one  flows  back  to  the 
condensation  module  and  the  second  to  the  plate-type  heat  exchanger.  Figure  2.2(a)  and 
2.2(b)  show,  respectively,  a  photo  of  the  micro-channel  condensation  module  and  another 
of  the  main  part  of  the  experimental  set-up. 


2.2.  Micro-Channel  Condensation  Module 

Figures  2.3(a)  and  2.3(b)  show  the  construction  of  micro-channel  condensation 
module.  The  module  consists  of  a  cover  plate,  housing,  condensation  copper  block,  and 
water  channels.  The  cover  plate,  which  seals  the  tops  of  the  square  micro-channels,  is 
made  from  transparent  polycarbonate  (Lexan)  to  facilitate  photographic  study  of  the 
condensing  flow.  The  top  of  the  oxygen-free  copper  block  is  2-cm  wide  by  29.9-cm  long, 
and  contains  ten  of  1  x  1  mm2  square  micro-channels.  Soldered  to  the  underside  of  the 
copper  block  are  three  of  3.8  x  3.8  mm2  brass  tubes,  which  carry  the  counter  flow  of 
cooling  water.  Sixteen  pairs  of  type-E  thermocouples,  which  are  embedded  in  the  copper 
block  beneath  the  micro-channels,  run  down  the  length  of  the  copper  block  at  19  mm 
intervals.  Using  the  assumption  of  one-dimensional  heat  conduction  between  the  two 
thermocouple  planes,  both  the  heat  flux  and  surface  temperature  at  the  base  of  the  micro¬ 
channels  may  be  determined  along  the  stream-wise  direction  (Kim  and  Mudawar,  20 1 0a). 
The  dimensions  of  the  test  section  are  provided  in  Table  2.2. 

The  copper  block  is  inserted  into  an  insulating  G-10  housing  that  features  FC-72 
inlet  and  outlet  ports,  micro-channel  inlet  and  outlet  plenums,  and  pressure  and 
temperature  measurement  ports.  The  G-10  housing  and  water  channels  are  covered  by 
insulating  layers  of  fiberglass. 


Table  2.2  Test  section  dimensions. 
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Figure  2.2  Photos  of  (a)  micro-channel  condensation  module  and  (b)  main  part  of  test  setup. 
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Figure  2.3  (a)  Construction  and  (b)  cross-sectional  view  of  micro-channel  condensation  module. 
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2.3.  Operating  Conditions  and  Measurement  Accuracy 

Experiments  were  performed  to  examine  flow  regimes  and  measure  pressure  drop 
and  heat  transfer  coefficients  corresponding  to  different  flow  rates  of  FC-72  and  cooling 
water.  The  test  matrix  for  the  study  consisted  of  24  tests  covering  six  FC-72  flow  rates 
(G  =  68,  1 18,  186,  248,  306,  and  367  kg/m2s)  and  four  water  flow  rates  (  mw  =  3,  4,  5,  and 
6  g/s). 

Since  the  focus  of  this  study  is  high-flux  condensation,  operating  conditions  were 
set  to  achieve  mostly  annular  flow  near  the  inlet  of  the  micro-channels.  Therefore,  the 
FC-72  was  introduced  into  the  condensation  module  in  pure  vapor  state  with  a  quality  of 
1.1 1  -  1.17  for  low  FC-72  flow  rates  (<7  =  68  -  186  kg/m2s)  and  1.08  -  1.10  for  high  FC- 
72  flow  rates  ( G  =  248  -  367  kg/m2s).  These  conditions  resulted  in  an  upstream  single¬ 
phase  superheated  vapor  region  1.1  -  3.1  cm  long  (4  -  10%  of  the  total  channel  length)  for 
the  low  FC-72  flow  rate  range,  and  2.7  -  5.8  cm  long  (9  -  19%  of  the  total  channel  length) 
for  the  high  FC-72  flow  rate  range.  The  inlet  plenum  pressure  ranged  from  1.040  -  1.324 
bar.  Operating  conditions  along  the  channel  length  were  as  follows:  saturation 
temperatures  of  Tsat  =  57.2  -  62.3°C,  qualities  of  x  =  0  -  1.17,  and  heat  fluxes  of  q"K  = 
0.43-3.21  W/cm2. 

There  are  several  techniques  for  measuring  liquid  film  thickness  in  two-phase 
applications  (Koskie  et  al.,  1989;  Lyu  and  Mudawar,  1991a,  1991b).  Unfortunately, 
small  channel  size  precludes  the  use  of  these  techniques  in  the  present  application, 
causing  a  reliance  on  photographic  methods  instead.  High-speed  video  imaging  played  a 
crucial  role  in  capturing  two-phase  condensation  flow  regimes.  Two  key  requirements 
for  capturing  the  complex  interfacial  features  in  the  micro-channels  with  high  resolution 
are  high  shutter  speed  and  high  magnification.  These  goals  were  achieved  with  a  Photron 
FASTCAM-Ultima  camera  capable  of  shutter  speeds  up  to  1/120,000  s,  which  was  used 
in  conjunction  with  an  assortment  of  Infinity  K-2  lenses. 

FC-72  pressure  was  measured  in  the  inlet  plenum  with  the  aide  of  an  absolute 
pressure  transducer,  while  a  differential  pressure  transducer  measured  pressure  drop 
between  the  inlet  and  outlet  plenums.  Temperatures  in  the  inlet  and  outlet  plenums, 
copper  block,  and  water  channels  were  measured  and  processed  by  an  HP-3852A  data 
acquisition  system.  All  measurements  were  made  after  all  the  system  pressures  and 
temperatures  reached  steady  state. 


Boiling  and  Two-Phase  Flow  Laboratory 


14 


By  using  five  thermocouple  ports  attached  on  top  of  the  cover  plate,  heat  loss 
through  the  cover  plate  is  estimated  at  less  than  2%  of  the  heat  input  through  the  base  of 
the  micro-channels,  which  is  calculated  using  the  assumption  of  one-dimensional  heat 
conduction  between  the  two  thermocouple  planes  (Kim  and  Mudawar,  2010a).  All 
pressure  transducers  used  in  this  study  featured  an  accuracy  of  0.05%.  The  rotameters 
had  a  flow  rate  accuracy  of  ±2%  of  measurement. 
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CHAPTER  3.  EXPERIMENTAL  RESULTS 

3.1.  Flow  Visualization  Results 

High-speed  video  imaging  and  photomicrographic  techniques  were  used  to 
capture  dominant  condensation  flow  regimes  along  the  parallel  micro-channels.  Figure 
3.1  depicts  representative  flow  regimes  in  three  adjacent  micro-channels  channels  at 
different  locations  from  the  micro-channel  inlet.  Five  distinct  flow  regimes  were 
identified.  The  smooth-annular  flow  regime  is  characterized  by  a  very  thin  and  fairly 
smooth  liquid  film  that  flows  along  the  channel  wall,  with  vapor  flowing  in  the  core  clear 
of  any  liquid  droplets.  The  film  is  barely  discernible  even  when  using  a  high 
magnification  lens.  This  regime  occurred  in  the  inlet  region  of  micro-channels  under 
conditions  of  high  inlet  quality.  The  wavy-annular  regime  features  a  liquid  film  that  is 
notably  thicker  that  that  of  the  smooth-annular  regime,  and  having  discernible  interfacial 
waves.  The  next  transition  regime  is  characterized  by  bridging  of  liquid  ligaments  across 
the  vapor  care.  The  slug  regime  features  elongated  cylindrical  bubbles  whose  length  is 
several  times  larger  than  the  width  of  the  channel.  The  bubbly  regime  features  spherical 
bubbles  with  a  diameter  approaching  the  of  the  micro-channel  width.  The  bubbly 
regimes  was  least  prevalent,  observed  only  at  the  lowest  tested  mass  velocity  of  G  =  68 
kg/m2s.  Notice  that  no  droplets  are  entrained  in  the  vapor  core  for  the  smooth-annular 
and  wavy-annular  flows.  This  is  in  sharp  contrast  with  annular  flow  associated  with  flow 
boiling  in  micro-channels.  As  observed  by  Qu  and  Mudawar  (2004),  liquid  droplets 
entrained  in  the  vapor  core  have  an  appreciable  influence  on  pressure  drop  and  heat 
transfer  in  micro-channel  boiling  flows.  This  fundamental  difference  between  micro- 
channel  flow  boiling  and  condensing  flows  and  its  implications  to  pressure  drop 
prediction  will  be  discussed  later. 

Figure  3.2(a)  shows  the  observed  FC-72  flow  regime  transition  data.  Increasing 
the  mass  velocity  is  shown  causing  the  smooth-annular  regime,  which  corresponds  to 
high  quality  values,  to  extend  further  downstream  towards  lower  quality  values,  while 
narrowing  the  slug  flow  in  the  low  quality  region.  Although  the  bubbly  regime  was 
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Figure  3.1  Representative  (a)  photographs  and  (b)  schematics  of  FC-72  condensation 
flow  regimes  in  three  adjacent  channels  for  G  =  68  kg/m2s. 
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(a) 


(b) 

Figure  3.2  Comparison  of  present  FC-72  condensation  flow  regime  data  with  those  of  (a) 
Wang  et  al.  (2002)  and  (b)  Triplett  et  al.  (1999)  and  Chung  and  Kawaji  (2004). 
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observed  in  the  present  study  for  G  =  68  kg/m2s,  no  bubbly  regime  data  are  shown  in  Fig. 
3.2(a).  This  is  due  to  the  fact  that  the  bubbly  regime  was  observed  for  negative  values  of 
quality,  indicating  bubbles  persisted  in  subcooled  liquid  flow  due  to  thermodynamic  non¬ 
equilibrium  effects.  Figure  3.2(a)  also  compares  the  present  FC-72  flow  regime 
boundary  data  with  flow  regime  boundary  lines  from  Wang  et  al.  (2002)  for  R134a 
condensing  flow  inside  a  horizontal  multi-port  condenser  featuring  rectangular  channels 
with  a  hydraulic  diameter  of  1 .46  mm.  Notice  that  both  the  smooth-annular  and  wavy- 
annular  regimes  identified  in  the  present  study  were  designated  simply  as  annular  by 
Wang  et  al..  The  boundary  between  the  transition  and  slug  regimes  of  the  present  study 
is  in  good  agreement  with  that  of  Wang  et  al.,  however,  the  boundary  between  the  wavy- 
annular  and  transition  regimes  corresponds  to  higher  quality  values  than  observed  by 
Wang  et  al.. 

Figure  3.2(b)  compares  the  present  FC-72  flow  regime  transition  data  with  those 
of  Triplett  et  al.  (1999)  and  Chung  and  Kawaji  (2004).  Triplett  et  al.’s  data  are  based  on 
air-water  adiabatic  flow  in  horizontal  semi-triangular  and  circular  micro-channels  with 
hydraulic  diameters  of  1.09  and  1.097  mm,  respectively.  Chung  and  Kawaji  conducted 
their  experiments  with  an  adiabatic  mixture  of  nitrogen  gas  and  water  inside  a  0.53-mm 
diameter  horizontal  circular  tube.  For  all  cases,  the  transition  line  to  slug  flow  is 
observed  when  the  superficial  gas  velocity  is  in  a  range  of  1  to  4  m/s.  The  transition  flow 
regime  of  this  study  was  designated  as  slug-annular  by  both  Triplett  et  al.  and  Chung  and 
Kawaji,  or  ‘serpentine-like  gas  core  in  churn  flow’  by  Chung  and  Kawaji.  The  wavy- 
annular  flow  regime  of  this  study  had  a  thicker  liquid  film  with  waves  having  larger 
amplitude  and  wavelength  than  those  of  Triplett  et  al.  and  Chung  and  Kawaji,  and  is 
similar  to  the  ‘gas  core  with  a  wavy  liquid  film’  regime  designated  as  churn  flow  by 
Chung  and  Kawaji.  Taking  into  consideration  the  differences  in  flow  regime  designation, 
the  present  flow  regime  boundaries  are  generally  consistent  with  those  from  the  two 
earlier  studies. 

Despite  the  fair  agreement  depicted  in  Figs.  3.2(a)  and  3.2(b),  it  is  important  to 
emphasize  the  fundamental  weaknesses  of  flow  regime  maps  that  utilize  dimensional 
plots.  The  first  is  that  the  use  of  G  versus  x  or  .7/  versus  Js  may  imply  that  such  pairs  of 
parameters  govern  all  flow  regimes,  which  is  obviously  not  the  case.  The  second 
weakness  is  the  fact  that  flow  regime  boundaries  are  dominated  by  physical  mechanisms 
and  corresponding  dimensionless  groups  that  vary  greatly  from  one  boundary  to  another. 


Boiling  and  Two-Phase  Flow  Laboratory 


19 


Figure  3.3  shows  an  alternative  presentation  of  flow  regime  boundaries  in  a 
dimensionless  plot  of  modified  Weber  number.  We  ,  versus  the  Martinelli  parameter,  X„. 
Using  the  assumption  that  the  inertia  of  the  vapor  phase  is  the  dominant  destructive  force 
acting  on  the  liquid  film,  while  surface  tension  and  liquid  viscous  forces  are  the 
stabilizing  forces,  Soliman  (1986)  derived  the  following  relations  for  modified  Weber 
number  by  balancing  the  destructive  and  stabilizing  forces. 


We  =  2.45 


Re 


0.64 


5«J3(l  +  1.09X“039)04 


for  Re,  <  1250, 
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where  the  Suratman  number,  Sus,  and  Martinelli  parameter  (based  on  turbulent  liquid- 
turbulent  vapor)  are  given  by 
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respectively  (Soliman,  1986).  Based  on  his  own  database,  Soliman  proposed  that  the 
flow  is  always  annular  for  We*  <  20  and  always  mist  for  We*  >  30.  While  a  few  of  the 
present  smooth-annular  flow  regime  data  fall  into  Soliman’s  annular-to-mist  transition 
region,  most  of  the  present  data  fall  in  Soliman’s  annular  regime  region.  Figure  3.3  also 
shows  boundaries  between  annular,  stratified-wavy,  and  plug  regimes  identified  by  Chen 
et  al.  (2006)  for  condensation  of  R134a  inside  horizontal  12  and  14  mm  diameter  micro- 
finned  tubes.  Overall,  there  is  general  agreement  with  the  present  boundary  lines 
between  the  wavy-annular  and  transition  regimes,  and  between  the  transition  and  slug 
regimes. 
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Figure  3.3  Comparison  of  present  FC-72  condensation  flow  regime  boundaries  with  those 

of  Soliman  (1986)  and  Chen  etal.  (2006). 
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3.2.  Two-Phase  Pressure  Drop  Results 


3.2.1.  Experimental  Results 

Figure  3.4  shows  the  variation  of  total  pressure  drop  measured  between  the  inlet 
and  outlet  plenums  of  the  condensation  module  with  water  mass  flow  rate  for  different 
FC-72  mass  velocities.  Increasing  the  mass  velocity  of  FC-72  increases  the  interfacial 
shear  stress  between  the  vapor  core  and  liquid  film  due  to  the  increasing  vapor  core 
velocity,  which  results  in  the  measured  increase  in  pressure  drop.  Figure  3.4  shows  the 
total  pressure  drop  decreases  slightly  with  increasing  water  mass  flow  rate.  This  trend 
can  be  explained  by  the  higher  water  flow  rate  increasing  the  wall  heat  flux,  which 
hastens  the  flow  deceleration  leading  to  diminished  interfacial  shear. 


3.2.2.  Pressure  Drop  Components 

For  a  superheated  inlet  condition,  the  total  pressure  drop  between  the  upstream 
and  downstream  plenums  can  be  obtained  from 

P*+AP*J+  (3-3> 

The  contraction  pressure  loss  and  expansion  recovery  at  the  inlet  and  outlet  of  the 
micro-channels,  respectively,  are  determined  from  relations  by  Collier  and  Thome  (1994), 
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For  single-phase  vapor  flow  at  the  plenum  inlet,  the  contraction  coefficient  Cc  due 
to  the  vena-contracta  is  obtained  from  a  relation  by  Geiger  (1964), 

1  — <7 
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2.08(l -®J+ 0.5371 ' 

For  single-phase  hydrodynamically  developing  vapor  flow  in  the  inlet  region  of 
the  micro-channels,  the  pressure  drop  can  be  expressed  as  a  function  of  the  apparent 
friction  factor. 
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Since  the  FC-72  is  introduced  into  the  condensation  module  in  superheated 
turbulent  state  with  a  Reynolds  number  range  of  6000  -  31,000,  the  local  pressure  is 
predicted  according  to  Zhi-qing  (1982)’s  analytical  solution  for  turbulent  developing 
vapor  flow  in  a  circular  tube, 

Lsp/Dh  =1.4039ReJf<Tl2?(l  +  0.1577<r  -0.1793<?+2  -0.0168<T3  +  0.0064<T4),  (3.7) 
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for  S*  =  1  (fully-developed  region).  (3.8b) 
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When  the  FC-72  completely  condenses  to  pure  liquid  in  the  downstream  region  of 
the  micro-channel,  the  pressure  drop  for  the  single-phase  liquid  region,  where  the  flow  is 
assumed  fully-developed,  can  be  determined  from  (Shah  and  London,  1978;  Incropera 
and  Dewitt,  2002) 
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for  Resp  /  <  2,000, 

(3.10a) 

fspJ  =  0.079Re;^5  for  2,000  <  Re.^  <  20,000, 

(3.10b) 

and 

LP,f  =0.046ReJj  for  Re^  >  20,000. 

(3.10c) 

For  the  two-phase  region,  the  pressure  drop  consists  of  accelerational  and 
frictional  components, 

U>Ip=*PIp,a+AP'P,f-  (3-H) 

For  a  constant  mass  velocity,  the  differential  form  of  accelerational  two-phase 
pressure  drop  can  be  expressed  as 


(3.12) 


r  dp N 

-G2  — 

vgx2  1  uf(\-xf 

[dz  j 

A  dz 

a  (l-«) 

where  the  void  fraction  can  be  obtained  from  Zivi  (1964)’s  correlation 
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1-x 


X 


V^3 


-1-1 


(3.13) 


For  the  homogeneous  equilibrium  model,  the  void  fraction  is  related  to 
thermodynamic  equilibrium  quality  by  the  relation 


a  = 


1  + 

1  \ 

l-x 

IV 

{  X  ) 

V 

-i-i 


(3.14) 


Substituting  the  above  relation  in  Eq.  (3.12)  while  neglecting  property  changes  yields  the 
following  widely  used  form  of  accelerational  two-phase  pressure  gradient  for  the 
homogeneous  equilibrium  model, 


(dP\  dx 

—  =Gvh — • 

dz  )A  '  dz 


(3.15) 


The  accelerational  two-phase  pressure  drop  can  be  determined  by  integrating  Eqs. 
(3.12)  or  (3.15)  numerically  along  the  streamwise  direction, 


A  P, 


fp,A 


(3.16) 


When  determining  the  total  pressure  drop  via  the  homogeneous  equilibrium 
models,  both  the  simple  and  general  forms  of  accelerational  pressure  drop,  Eqs.  (3.15) 
and  (3.12),  respectively,  will  be  examined.  In  case  of  separated  flow  models,  only  the 
general  form  of  accelerational  pressure  gradient,  Eq.  (3.12),  will  be  used  to  determine  the 
accelerational  pressure  drop. 


The  frictional  two-phase  pressure  drop  can  be  determined  by  integrating  the 
relation  for  frictional  pressure  gradient  corresponding  to  the  chosen  model  or  correlation 
as  indicated  in  Tables  3.1  -  3.3, 


A P, 


tp.F 


dz . 

F 


(3.17) 


In  the  two-phase  condensing  region,  the  local  saturation  temperature  of  FC-72  is 
obtained  from  the  local  saturation  pressure,  and  the  thermophysical  properties  for  liquid 
and  vapor  are  based  on  local  saturation  pressure. 
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3.2.3.  Comparison  with  Different  Pressure  Drop  Models  and  Correlations 

To  determine  the  two-phase  frictional  pressure  drop,  both  homogeneous 
equilibrium  and  separated  flow  models  are  considered. 

Homogeneous  Equilibrium  Model  (HEM) 

Both  constant  friction  factor  (Bottomley,  1936;  Lewis  and  Robertson,  1940; 
Benjamin  and  Miller,  1942;  Markson  et  al.,  1942;  Allen,  1951)  and  two-phase  mixture 
viscosity  (McAdams  et  al.,  1942;  Akers  et  al.,  1958;  Cicchitti  et  al.,  1960;  Owens,  1961; 
Dukler  et  al.,  1964;  Beattie  and  Whalley,  1982;  Lin  et  al.,  1991)  methods  are  used  to 
predict  the  two-phase  frictional  pressure  gradient  via  the  homogeneous  equilibrium 
model.  In  order  to  check  the  sensitivity  of  different  accelerational  pressure  drop  relations 
to  total  pressure  drop,  both  the  simple  and  general  forms  of  accelerational  pressure  drop, 
Eqs.  (3.15)  and  (3.12),  respectively,  are  compared.  Figures  3.5(a)  and  3.5(b)  compare  the 
present  FC-72  pressure  drop  data  to  predictions  based  on  two  different  constant  friction 
factors,  f,p  =  0.003  and  0.005,  and  seven  different  two-phase  viscosity  models.  The  total 
pressure  drop  in  Fig.  3.5(a)  is  based  entirely  on  the  homogeneous  equilibrium  model  (i.e., 
using  Eq.  (3.15)  for  the  accelerational  pressure  gradient),  while  the  accelerational 
pressure  gradient  in  Fig.  3.5(b)  is  based  on  Eq.  (3.12)  and  Zivi’s  relation  for  void  fraction. 
Comparing  Figs.  3.5(a)  and  3.5(b)  shows  the  former  approach  only  slightly  overpredicts 
the  data  compared  to  the  latter.  Therefore,  the  differences  in  predictions  of  total  pressure 
drop  in  Figs.  3.5(a)  and  3.5(b)  are  due  mostly  to  the  different  frictional  pressure  drop 
models  used. 

Given  the  slightly  better  predictions  in  Fig.  3.5(b)  compared  to  Fig.  3.5(a),  all 
subsequent  accelerational  pressure  drop  calculations  in  Figs.  3.6,  3.7(a)  and  3.7(b)  will  be 
based  on  the  general  form  of  accelerational  pressure  drop,  Eq.  (3.12),  using  Zivi’s  void 
fraction  relation. 

Table  3.1  provides  details  of  the  constant  two-phase  friction  factor  and  two-phase 
viscosity  methods  as  well  as  their  accuracy  in  predicting  the  present  data.  The  accuracy 
of  individual  models  is  measured  by  mean  absolute  error,  which  is  defined  as 

i  UP  —  A P  I 

MAE  =  —Y x  1 00% .  (3.18) 

M  ^  A Pexp 

In  general,  the  homogeneous  equilibrium  model  is  better  suited  to  dispersed  two- 
phase  flows  such  as  bubbly  flow.  The  fundamental  premises  of  this  model  fall  apart  for 
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Table  3.1  Homogeneous  equilibrium  model  for  two-phase  frictional  pressure  gradient, 
and  corresponding  MAE  in  predicting  present  total  pressure  drop  data. 


=-7T^Pul  =  -TT^°tGiX  +  x— 

dz  )F  D.  '  D,  u, 

\  / 1  h  h  } 

Constant  friction  factor  method 

Author(s) 

Applications 

fp 

MAE  [%] 

Lewis  and  Robertson  ( 1 940) 
Markson  e/ r;/.  (1942) 

High  pressure  steam-water  boilers 

0.005 

31.72 

Botlomley  (1936) 

Benjamin  and  Miller  ( 1 942) 
Allen  (1951) 

Low  pressure  flashing  steam-water  flows 

0.003 

56.77 

Two-phase  mixture  viscosity  method 

fp  Re,,,  =  24(l  -1.3553)?  + 1.9467/?2  -1.7012/?'  +  0.9564/?4  -0.2537/95)  for  Re, „  <2,000 

4  =  0.079  Re,; 25  for  2,000  <  Re,„  <  20,000 

4  =  0.046  Re,;2  for  Re,„  >  20,000  ,  where  Re,„  =  GDh/^ 

Author(s) 

Equation 

MAE  [%] 

McAdams  et  at.  (1942) 

1  x  +  l-x 

1  "7 

P.p  F  P, 

1  /  .2.1 

Akers  et  al.  (1958) 

P'p  = 

Pf 

34.98 

1  v 

(1-4  +  x  -£- 

KJ 

Cicchitti  et  al.  (1960) 

A- =  x  A  +<l 

85.16 

Owens  (1961) 

11 

306.0 

Duklere/a/.  (1964) 

_«V*,+0--*)y/A/ 

XUg  +  (l-x)L)f 

15.92 

Beattie  and  Whalley  (1982) 

P,„  =  (Ofig+{\-o))(\  +  2.5co)fjf 

xv 

(O  = - - — 

Uf+XUK 

17.34 

Lin  et  al.  (1991) 

P/Ps 

Pg  +x'4(Pf~Pg) 

17.05 
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Figure  3.5  Comparison  of  present  FC-72  pressure  drop  data  with  predictions  of 
(a)  homogeneous  equilibrium  model,  and  (b)  homogeneous  equilibrium  model 
with  accelerational  pressure  drop  based  on  separated  flow  model  and 
Zivi  (1964)’s  relation  for  void  fraction. 
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regimes  involving  separation  between  the  liquid  and  vapor  phases  or  large  velocity 
differences  between  the  two  phases.  Annular  flow  is  an  example  of  a  regime  where  the 
use  of  the  homogenous  equilibrium  model  is  generally  avoided.  In  spite  of  the  fact  that 
annular  flow  was  dominant  throughout  the  present  measurement,  Figs.  3.5(a)  and  3.5(b) 
prove  that,  when  used  in  conjunction  with  most  two-phase  viscosity  models,  excepting 
those  of  Cicchitti  et  al.  (1960)  and  Owens  (1961),  the  homogeneous  equilibrium  model 
provides  fairly  accurate  predictions  of  the  present  experimental  data.  Most  notably,  the 
viscosity  model  by  Dukler  et  al.  (1964),  with  a  MAE  of  15.92%,  provides  the  best 
predictions  among  all  viscosity  relations.  In  fact,  this  model  provides  better  accuracy 
than  all  separated  flow  correlations  as  discussed  in  the  next  section. 

Figure  3.6  shows  the  contributions  of  individual  pressure  drop  components  to 
total  pressure  drop  using  the  mixture  viscosity  model  by  Dukler  et  al..  For  a  condensing 
flow,  the  accelerational  two-phase  pressure  drop  term  is  negative  because  of  the 
deceleration  along  the  stream-wise  direction.  In  fact,  condensation  reduces  two-phase 
pressure  drop,  as  opposed  to  flow  boiling,  where  two-phase  pressure  drop  is  increased  by 
the  stream-wise  acceleration,  and  the  accelerational  pressure  drop  is  positive.  Figure  3.6 
shows  that  the  magnitude  of  AP,0,  is  dictated  mostly  by  two-phase  friction  and 
acceleration.  Increasing  the  mass  velocity  decreases  the  magnitude  of  A Plp  A  /  A Plol 
because  of  the  decreasing  amount  of  vapor  being  condensed.  The  contributions  of  A Pspf 
and  A Pt,  to  AP„„  are  minimal,  while  those  of  A Pc  and  A Pspd  are  more  significant  because 
of  the  pure  vapor  state  at  the  micro-channel  inlet. 

Separated  Flow  Models  (SFMs) 

Tables  3.2  and  3.3  provide  a  select  summary  of  separated  flow  frictional  pressure 
drop  correlations  that  have  been  recommended  for  macro-channels  (Lockhart  and 
Martinelli,  1949;  Friedel,  1979;  Chisholm,  1983;  Muller-Steinhagen  and  Heck,  1986; 
Jung  and  Radermacher,  1989;  Wang  et  al.,  1997;  Cavallini  et  al.,  2002)  and  mini/micro¬ 
channels  (Mishima  and  Hibiki,  1996;  Tran  et  al.,  2000;  Chen  et  al.,  2001;  Lee  and  Lee, 
2001;  Zhang  and  Webb,  2001;  Yu  et  al.,  2002;  Lee  and  Mudawar,  2005;  Nino  et  al., 
2005;  Hwang  and  Kim,  2006),  respectively.  These  include  a  mix  of  correlations  for  flow 
boiling  and  condensation,  adiabatic,  flow  boiling  and  condensing  flows,  and  for  isolated 
tubes  and  multi-channel  heat  sinks.  Notice  that  the  correlation  of  Jung  and  Radermacher 
( 1 989)  is  based  on  the  assumption  of  A Plp  «  A P,  because  the  accelerational  component 
was  found  to  be  small  in  their  flow  boiling  experiments. 
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Figure  3.6  Contributions  of  individual  components  of  pressure  drop  to  total  pressure  drop 
for  different  mass  velocities  and  mw  =  3g/s  predicted  with  the  frictional  pressure  gradient 
determined  by  the  homogenous  equilibrium  model  using  the  viscosity  model  of 
Dukler  etal.  (1964),  and  the  accelerational  pressure  gradient  determined  by 
the  separated  flow  model  using  Zivi  (1964)’s  relation  for  void  fraction. 
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Table  3.2  Two-phase  frictional  pressure  gradient  correlations  for  macro-channels  based  on  the 
separated  flow  model,  and  corresponding  MAE  in  predicting  present  total  pressure  drop  data. 


Author(s) 

Equation 

Remarks 

MAE 

[%] 

Lockhart  and 
Martinelli 
(1949) 

\dz)F  \dz)f'  '  X  X 

C„  =  5,  C„.  =  10,  C„=12,  C„  =  20 

(dp/dz)g 

Dh=  1.49-25.83  mm 
Adiabatic 

Fluids:  water,  oils, 
hydrocarbons 

80.26 

Friedel 

(1979) 

[f]  =ffl 

V  dz  )F  v  dz  )fo 

(  V’-9 

+3.24x°'78  (l  -  x)°214 1  — 

Fr=  °2  We=G2D* 

S  L>h  Ph  ’  °P„ 

-  x)2  +  x 2 

1  >0.19 

A 

U/J 

•  =  — 
XU 

w«¥4.)+ 

y/JUJ 

\°-7 

1-&-  Fr,-;M5We?035 

Pf) 

1 

,+0 -■*)«/ 

D  >  4  mm 

Fluids:  air/water,  air/oil, 
R12  (25000  data  points) 

173.4 

Chisholm 

(1983) 

(f  1  =  (f ),  *  •  * = 1 +<r'  -1)[s/'” (1  “  9'" + «■■’] 

(dP/dz) 

T2  =  - - ,  for  B  values  refer  to  Chisholm  (1983) 

(dPldz)lo 

Flow  boiling 

Fluids:  air/water,  steam 

114.1 

Muller- 
Steinhagen 
and  Heck 
(1986) 

( iip)  _  \( iip)  +9 

V  dz  Jr  (l  dz  jfo  \dz)g0 

-(fi; 

]  ,,  u/3  f  rtf’ 3  5 

T  > 

D  =  4  -  392  mm 

Fluids:  air/water,  water, 
hydrocarbons, 
refrigerants 
(9300  data  points) 

54.88 

Jung  and 
Radermacher 
(1989) 

(—}  =  (— 1  f  f  =12.: 
1  dz  )F  Uz;„,  '  fo 

>0.9/  \0'5 

'•■fel  P?)  (?) 

82  AT,;1 47  (1 

-*)“ 

D  =  9. 1  mm 

Annular  flow  boiling 
Fluids:  pure  and 
mixtures  ofR22.  R114. 
R12.  R152a 

288.3 

Wang  et  al. 
(1997) 

For  G  >200  kg / m2 s ,  1  —  1  =(  —  1 

V  dz  )F  \  dz  )g 

</>;  =1  +  9.4  A062  +  0.564 X245 

For  G  <  200 kg/m2s ,  \ — )  =  (  — ]  )) 

V  dz  )F  \dz  ), 

V'V 

C  =  4.566  xlO-6*0128  Ref*  —  — 

V  A 

:,*;  =  i+£+-L, 

n  A-  A- 

\5-> 

D  =  6.5  mm 

Adiabatic 

Fluids:  R22.  R134a, 
R407C 

86.19 

Cavallini 
et  al.  (2002) 

f-1  =(-)  HL.We  -< 

\  dz  )F  l  dz  )/o  *°  , 

For/;* >2.5.  fl={ l-xf+x 

+1.262x0'6‘',8f — 1 

uf) 

For/;* <2.5, 

^,/* 
0..O-  * 

\/  r 

.2  ^  i_ 

Fr  JU 

0.3278  ✓ 

A 

U/J 

C.r  1  P.A 

A:  \(Pf  ~  Pg)sDh 

+... 

t,  J 

-1.181  x  x  3.477 

i-Ar.  1  IFc'-01458 

"/  J  8" 

D  >  3.1  mm 

Condensation 

Fluids:  R22.  R134a, 

R125,  R32,  R236ea, 
R407c,  R410A 

597.8 
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Table  3.3  Two-phase  frictional  pressure  gradient  correlations  for  mini/micro-channels  based  on 
the  separated  flow  model,  and  corresponding  MAE  in  predicting  present  total  pressure  drop  data. 


Author(s) 


Equation 


Remarks 


MAE 

[%] 


Mishima  and 
Hibiki 
(1996) 


dP)  ,,  ,2  ,  C  1 
dz),  '  '  X  X 


C  =  21[l-exp(-0.319Z)„)];D„  [mm] 


D  =  1.05-4. OS  mm 
Adiabatic 
Fluid:  air/water 


27.15 


Tran  et  al. 
(2000) 


JdP)  ,2  M 

dz\o  J«(P/ 


^  =  i+(4.3r2-i)[^-(i-xp+x'75] 


Dh  =  2.40.  2.46,  2.92  mm 
Flow  boiling 
Fluids:  R134a.  R12, 

R113 


276.1 


Lee  and  Lee 
(2001) 


dP)  JdP 
dz  )F  (  dz 


f  ^  =  i+£+_L,  v,=nj!  x  = 

f  '  XX-  a 


Pi 


P,°D„ 


Dh  =  0.78-6.67  mm 
Adiabatic 
Fluid:  air/water 


C„.  =6.833x10  /1 

1-2 


,-S  J -1.317, .,0.719 


Re,„ 


C„,=  3.627  Re  „ 


C„  =6.185  xl0~2  Re/726 


,  C„  =0.048  Re,, 


36.26 


Zhang  and 
Webb 
(2001) 


(dP 

1  _ 

(dP\ 

lit. 

[liz ) 

(>- 

,y)~  +  ^ 

fdP 

1  _ 

fdP') 

Ur, 

l  dz  J 

D  =  0.96-6. 20mm 
Adiabatic 
Multi-channel 
Fluids:  R134a,  R22, 
R404a 


522.9 


Chen  et  al. 
(2001) 


Q 

fo.Friedel 

0.0333  Reu; 45 

For  Bo  <2.5,  Sl  =  - 


D  =  1.02-9  mm 
Adiabatic 
Fluids:  air/water, 
R410A.  ammonia 


Re1]1"' (l  + 0.4c  j 


For  Bo  >  2.5 ,  £2  = 


We!„ 


(2.5  +  0.06  Bo) 


■  Bo  =  g(p,  -  a,) 


(Djf 


47.65 


Y  u  et  al. 
(2002) 


(dP)  _ 

(  dz  Jf. 

[dz], 

xi 


X.  =18.65 


> 

°.5  ,  ,  . 

J 

l  ^  Ji 

Re" 


D  =  2.98  mm 
Flow  boiling 
Fluids:  water,  ethylene 
glycol,  aqueous 
mixtures  of  ethylene 
glycol 


51.25 


Nino  et  al. 
(2005) 


(in 

l  dz)p 

( dz  K- 

it  ,  We,  = 


aP* 


J  =  exp(-0.046A„„, )  +  0.22[exp(-0.002Yu„„)  -  exp(-7A,„„ )] 


Dh=  1.02,  1.54  mm 
Adiabatic 
Multi-channel 
Fluids:  R410A,  R134a, 
air/water 


X . = 


V  \ 

x  +  1 

iV 

0.9  " 

y  — 

' 

-  0,,5 

V| 

"  We'/ 

it  - 

Uw 

brJ 

V  v« , 

34.24 


Lee  and 
Mudawar 
(2005) 


(in 

)  JdP) 

(  dz  . 

l  (dz) 

C,  = 

2-16Re/ 

(in 

|  JdP) 

(  dz  . 

l  (dz) 

Vf  >  Vf 


=  2. 1 6  Re  fo0047  We°/°  ,  C„  =  1 .45  Re  We'J  ,  Wefo  = 


^;='+T7 


C  1 

X  X~ 


GJ\ 

ap, 


Dh  =  0.349  mm 
Flow  boiling 
Multi-channel 
Fluids:  R134a,  water 


57.55 


Hwang  and 
Kim  (2006) 


,  /)  =  l  +  £.  +  _L ,  C  =  0.227  Re"f 2  A"0  33 A';082 
XX1 


l;,/} 


conf 


D  =  0.244-0.792  mm 

Adiabatic 

Fluid:  R134a 


53.95 
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Figure  3.7  Comparison  of  present  FC-72  pressure  drop  data  with  predictions  of  separated 
flow  correlations  recommended  for  (a)  macro-channels  and  (b)  mini/micro-channels. 
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Figure  3.7(a)  shows  that  all  seven  macro-channel  correlations  overpredict  the 
present  FC-72  condensation  data.  Notably,  the  total  pressure  drop  is  highly  overpredicted 
by  the  correlations  of  Friedel  (1979),  Chisholm  (1983),  Jung  and  Radermacher  (1989), 
and  Cavallini  et  al.  (2002),  with  corresponding  MAE  values  of  173.42%,  114.06%, 
288.27%,  and  597.75%,  respectively. 

Figure  3.7(b)  compares  the  present  FC-72  micro-channel  condensation  pressure 
drop  data  with  predictions  using  separated  flow  correlations  for  frictional  pressure  drop  in 
mini/micro-channels.  The  correlations  of  Mishima  and  Flibiki  (1996),  Lee  and  Lee 
(2001),  and  Nino  et  al.  (2005),  which  are  all  based  on  adiabatic  air-water  flow,  provide 
fair  predictions  of  the  present  data,  with  MAE  values  of  27.15%,  36.26%,  and  34.24%, 
respectively.  The  correlations  of  Lee  and  Mudawar  (2005a)  and  Tran  et  al.  (2000), 
which  are  based  on  micro-channel  flow  boiling,  overpredict  the  present  data  with  a  MAE 
of  57.55%  and  276.1%,  respectively.  The  large  MAE  of  522.86%  for  Zhang  and  Webb 
(2001)’s  correlation  may  be  attributed  to  the  relatively  low  reduced  pressure  range  of  the 
present  FC-72  data,  Pr  =  0.05-0.07,  compared  to  the  recommended  range  of  Pr  >  0.2  for 
the  original  correlation. 

By  closely  examining  the  predictive  accuracy  of  the  separated  mini/micro¬ 
channels  flow  correlations,  it  can  be  concluded  that,  even  on  a  local  basis,  condensation 
behavior  is  fundamentally  different  from  that  of  flow  boiling  and  perhaps  closer  to  that  of 
adiabatic  two-phase  flow.  One  obvious  difference  between  flow  boiling  on  one  hand, 
and  condensing  and  adiabatic  two-phase  flows  on  the  other,  is  the  relative  abundance  of 
entrained  droplet  in  the  former  and  absence  from  the  latter  two,  as  schematically  shown 
in  Fig.  3.8(a).  This  phenomenon  is  especially  important  for  mini/micro-channel  flows, 
where  the  annular  regime  is  dominant  regardless  of  heating  conditions.  As  indicated  by 
Qu  and  Mudawar  (2004),  the  transition  to  the  annular  regime  in  flow  boiling  occurs  far 
upstream  in  a  micro-channel,  but  with  an  abundance  of  entrained  droplets  that  are  formed 
by  shattering  of  liquid  from  the  micro-channel’s  upstream.  However,  the  flow 
visualization  results  of  the  present  study  show  no  droplets  are  formed  in  the  annular 
regions  of  the  micro-channel.  This  conclusion  highlights  the  need  for  new  models  and 
correlations  that  can  accurately  predict  condensation  behavior  in  mini/micro-channels, 
especially  when  implemented  in  a  multi-channel  configuration. 


Boiling  and  Two-Phase  Flow  Laboratory 


34 


Vapor  Control 

core  volume 


Liquid  film 


Annular  Condensation 


Liquid  film 


Low  surface  tension  High  surface  tension 


Figure  3.8  (a)  Fundamental  differences  between  annular  condensation  and  annular  flow 
boiling  in  micro-channels,  (b)  Effects  of  surface  tension  on  liquid  distribution 

along  channel  perimeter. 
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3.3.  Condensation  Heat  Transfer  Results 


3.3.1.  Heat  Transfer  Data  Reduction 


Efforts  were  made  to  achieve  slightly  superheated  FC-72  conditions  at  the  inlet  to 
the  condensation  module.  This  helped  maintain  annular  flow  in  the  inlet  region  of  the 
micro-channels  and  achieve  high  heat  transfer  coefficients  as  well  as  a  broad  range  of 
operating  conditions. 


Figure  3.9  shows  a  unit  cell  illustrating  key  parameters  that  are  used  in 
determining  the  heat  transfer  performance  of  the  condensation  micro-channels.  Unlike 
most  boiling  experiments,  where  the  heat  flux  is  measured  from  electrical  power  input, 
heat  flux  in  this  study  is  determined  from  temperature  measurements  along  two  parallel 
planes  in  the  copper  block  beneath  the  micro-channels.  Sixteen  pairs  of  type-E 
thermocouples  are  embedded  at  19-mm  intervals  along  the  length  of  the  copper  block.  At 
each  axial  thermocouple  location,  both  the  heat  flux  along  the  base  of  the  micro-channels 
and  the  base  temperature  are  determined  from  temperature  measurements  by  the 
thermocouple  pair  using  the  assumption  of  one-dimensional  heat  conduction  between  the 
two  thermocouple  planes. 


V  base 


K(T,  -T>) 


H„ 


(3.19) 


The  temperature  corresponding  to  the  plane  in  line  with  the  bottom  of  the  micro¬ 
channels  can  be  determined  from 


Tw,b  =1  + 


Vl.se  H, 


(3.20) 


The  local  rates  of  sensible  heat  loss  in  the  superheated  vapor  region  and  latent 
heat  loss  in  the  two-phase  condensing  region  from  FC-72  can  be  calculated  from  the 
following  energy  balances,  respectively, 


jfCpjA  rf=  ql„  {Wcl+W,)Az,  (3.21  a) 

and  ^ha& x=ql,(WA+W,)bz.  (3.21b) 

The  onset  of  two-phase  condensation  occurs  at  a  distance  from  the  inlet  of  the 
micro-channel  where  the  fluid  temperature  reaches  the  saturation  temperature 
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corresponding  to  local  pressure.  The  local  quality  in  the  superheated  region  can  be 
determined  from 


x 


,g  {Tt  _ 


'i>,g 


‘fg 


(3.22) 


Figure  3.10  shows  a  sample  of  top  and  bottom  thermocouple  temperatures 
measured  along  the  copper  block  for  three  FC-72  mass  velocities.  The  corresponding 
temperature  distribution  along  the  two  thermocouple  planes  is  fitted  by  a  third-order  or 
fourth-order  polynomial  function  of  axial  distance.  These  distributions  are  then  used  to 
determine  the  corresponding  variations  of  base  heat  flux  and  thermodynamic  equilibrium 
quality.  Figure  3.10  shows  that  the  base  heat  flux  increases  with  increasing  mass  velocity 
due  to  the  increasing  temperature  difference  between  the  two  thermocouple  planes. 


Inlet  quality  was  maintained  at  1.11  -  1.17  for  low  FC-72  mass  velocities  (G  =  68 
-  186  kg/m2s)  and  1.08  -  1.10  for  high  mass  velocities  (G  =  248  -  367  kg/m2s),  which 
correspond  to  a  single-phase  superheated  vapor  region  1.1  -  3.1  cm  long  (4  -  10%  of  the 
channel  length)  and  2.7  -  5.8  cm  long  (9  -  19%  of  the  channel  length),  respectively.  For 
the  two-phase  condensation  region  (x  <  1),  the  present  study  covered  the  following 
operating  conditions:  FC-72  mass  velocities  of  G  =  68  -  367  kg/m2s,  FC-72  saturation 
temperatures  of  Tsat  =  57.2  -  62.3°C,  FC-72  qualities  of  x  =  1  -  0,  water  mass  flow  rates 
of  mw  =3-6  g/s,  and  heat  fluxes  of  q"v  =  0.43  -  3.21  W/cm2. 


3.3.2.  Experimental  Results 

Kim  and  Mudawar  (2010a,  2010b)  showed  that  the  one-dimensional  fin  analysis 
method  is  both  convenient  and  accurate  in  evaluating  the  local  heat  transfer  coefficient 
for  micro-channels  with  rectangular,  inverse  trapezoidal,  triangular,  trapezoidal, 
diamond-shaped,  and  circular  cross  sections. 

In  case  of  a  condenser  having  rectangular  micro-channels  with  three-sided  cooling 
walls  (/.<?.,  assuming  a  perfectly  insulating  top  cover  plate)  as  illustrated  in  Fig.  3.9, 
applying  the  fin  analysis  method  yields 

?:(2g„+yJ  _  qlM+Wj  .  . 

(T,J,-T,PnHcl1+W,t)  (T^-TjlnH^+W^)'  l~~ 

where  qbase  is  the  base  heat  flux,  and  the  fin  efficiency  and  fin  parameter  are  defined, 
respectively,  as  (Incropera  and  Dewitt,  2002) 
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Figure  3.10  Variations  of  copper  block  temperatures,  base  heat  flux,  and  thermodynamic 
equilibrium  quality  along  stream-wise  direction  for  different  FC-72  mass  velocities 

with  mw  =  3  g/s. 
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n=- 


tanh  (mHch ) 


mH 


ch 


(3.24a) 


and  in—  .  (3.24b) 

V.,ws 

Figure  3.11  shows  the  variation  of  the  experimentally-determined  local  FC-72 
condensation  heat  transfer  coefficient  with  thermodynamic  equilibrium  quality  in  the 
saturated  region  for  different  water  mass  flow  rates  and  FC-72  mass  velocities  in  the 
range  of  68  -  367  kg/m2s.  Very  high  heat  transfer  coefficient  values  are  achieved  near  x 
=  1,  where  annular  flow  regime  is  initiated  and  the  liquid  film  is  very  thin.  The 
condensation  coefficient  decreases  along  the  channel  as  the  film  gradually  thickens 
and/or  the  annular  regime  is  replaced  by  transition  or  slug  flow  regimes.  Overall,  there  is 
a  far  stronger  sensitivity  to  FC-72  mass  velocity  than  to  water  flow  rate. 


Figure  3.12  shows  the  variation  of  local  FC-72  condensation  heat  transfer 
coefficient  with  thermodynamic  equilibrium  quality  in  the  saturated  region  for  different 
FC-72  mass  velocities  and  two  water  flow  rates.  Increasing  the  mass  velocity  of  FC-72 
extends  the  annular  regime  farther  downstream  towards  lower  quality  values.  Increasing 
the  mass  velocity  increases  the  condensation  heat  transfer  coefficient  by  thinning  the 
liquid  film  due  to  the  increased  interfacial  vapor  shear  stress.  Notice  that  the  slope  of 
heat  transfer  coefficient  plot  in  Fig.  3.12  decreases  appreciably  in  range  of  x  =  0.4  -  0.7 
for  G  =  248  kg/m2s  and  a-  =  0.6  -  0.8  for  G  =  306  kg/m2s.  As  presented  in  Chapter  3.1, 
these  quality  ranges  correspond  mostly  to  the  wavy-annular  regime.  It  can  therefore  be 
inferred  that  at  high  mass  velocities,  interfacial  waves  enhance  the  condensation  heat 
transfer  by  slowing  the  trend  of  declining  heat  transfer  coefficient  along  the  channel 
caused  by  the  film  thinning.  Notice  that  the  flattening  in  the  variation  of  h,,,  versus  x  is 
nonexistent  for  the  lower  mass  velocities.  This  can  be  explained  by  the  narrowing  of  the 
extent  of  the  wavy  annular  regime  at  low  mass  velocities. 


Figure  3.13  shows  the  variation  of  the  average  condensation  heat  transfer 
coefficient  with  water  mass  flow  rate  for  different  FC-72  mass  velocities.  It  should  be 
emphasized  that  all  of  the  average  condensation  heat  transfer  coefficient  data  in  this 
figure,  as  well  as  Figs.  3.14-3.16  and  3.20,  correspond  to  values  averaged  over  only  the 
two-phase  condensing  region  (corresponding  to  the  range  x  =  0  -  1),  and  exclude  the 
upstream  superheated  vapor  region  and  downstream  single-phase  liquid  region  of  the 
micro-channel.  In  Fig.  3.13,  the  condensation  heat  transfer  coefficient  increases  with 
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Figure  3.12  Variation  of  experimentally-determined  local  FC-72  condensation  heat 
transfer  coefficient  with  thermodynamic  equilibrium  quality  for  different  FC-72 
mass  velocities  with  (a)  mw  =  3  g/s,  and  (b)  mw  =  6  g/s. 
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increasing  FC-72  mass  velocity  because  of  the  increased  interfacial  shear  stress  and 
resulting  liquid  film  thinning.  On  the  other  hand,  increasing  the  water  flow  rate  increases 
the  flow  rate  of  FC-72  and  film  thickness,  casing  a  reduction  in  the  heat  transfer 
coefficient. 
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Figure  3.13  Variation  of  experimentally-determined  average  FC-72  condensation  heat 
transfer  coefficient  with  coolant  mass  flow  rate  for  different  FC-72  mass  velocities. 
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3.3.3.  Comparison  with  Previous  Heat  Transfer  Correlations 


Table  3.4  provides  a  summary  of  select  condensation  heat  transfer  correlations  for 
the  annular  flow  regime.  The  predictive  accuracy  of  the  correlations  is  measured  by 
mean  absolute  error,  which  is  defined  as 


MAE  =  —  yl/l'p’p-  /'Ip’exp  x  100%.  (3.25) 

M  /?tP-e*P 

The  correlations  in  Table  3.4  can  be  classified  into  the  following  categories:  two- 
phase  multiplier-based  correlations  (Akers  et  al.,  1958;  Shah,  1979;  Dobson  and  Chato, 
1998),  and  boundary  layer-based  correlations  (Cavallini  and  Zecchin,  1974;  Moser  et  al., 
1998;  Wang  et  al.,  2002;  Koyama  et  al.,  2003;  Huang  et  al.,  2010).  Notice  that  the 
correlations  of  Koyama  et  al.  and  Huang  et  al.  are  based  on  Haraguchi  et  al.' s  (1994b) 
boundary  layer-based  correlation.  Since  the  correlations  in  Table  3.4  are  intended  for 
uniform  circumferential  cooling  in  circular  tubes  or  rectangular  channels  with  four-sided 
cooling,  a  multiplier  must  be  adopted  when  applying  these  correlations  to  condensation  in 
rectangular  micro-channels  with  three-sided  wall  cooling  (i.e.,  with  a  perfectly  insulating 
top  cover  plate).  Following  a  technique  adopted  by  Qu  and  Mudawar  (2003b)  and  Lee 
and  Mudawar  (2005)  to  correct  for  three-sided  heating  in  flow  boiling,  the  condensation 
heat  transfer  coefficient  for  three-sided  cooling  is  related  to  that  obtained  from  a 
correlation  for  uniform  circumferential  cooling  by  the  relation 


K  = 


r  Nu^ 

v  NlU  j 


ltp,cir  ’ 


(3.26) 


where  h,pxir  is  the  heat  transfer  correlation  based  on  uniform  circumferential  cooling,  and 
Nu 3  and  Nu4  are  Nusselt  numbers  for  thermally  developed  laminar  flow  with  three-sided 
and  four-sided  heat  transfer,  respectively,  (Shah  and  London,  1978) 

Nu3  =  8.235(l - 1.833/?  +  3.767/?2  -5.814/?3  +  5.361/?4  -2.0/?5),  (3.27a) 

and  Nu4  =  8.235(l  —  2.042/?  +  3.085/?2  -2A11  ^  +1.058/?4  -0.1 86/?5).  (3.27b) 

Figure  3.14  compares  the  present  FC-72  average  condensation  heat  transfer 
coefficient  data,  which  correspond  mostly  to  the  smooth-annular,  wavy-annular,  and 
transition  flow  regimes,  to  predictions  of  five  correlations  for  annular  flow  condensation 
in  macro-channels.  Four  of  the  correlations  show  good  to  fair  predictions  of  the  present 
data,  while  the  correlation  by  Akers  et  al.  (1958)  highly  over-predicts  the  present  data. 
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Table  3.4  Condensation  heat  transfer  correlations  for  annular  flow  regime  with 
corresponding  MAE  in  predicting  present  experimental  data. 


Author(s) 


Equation 


Remarks 


MAE 

_[%] 


Akers  et  al. 
(1958) 


KDb  _  (  Nui 


0.0265  Re"  8Pr'  5  for  Re  >50000 
Nu.  1  "  ' 


\Dh  _  f  Nu 3 


Fluids:  R12,  propane 
All  flow  regimes 
60  <  Nu  Pr^',!  <  400 


Re„„  =  G 


5.03 Re1 ,3Prl  5  for  Re  <50000 
Nu  ,  '  "  ' 


A 


■ 

\ 

(1  -x)  +  x| 

pf_ 

UJ 

A 


160.0 


Cavallini 

and 

Zecchin 

(1974) 


jp_  _| 
k, 


/  \ 

“ 

„  A 

0.5 

1  N‘h 

0.05Re"8PrOJ3 

'  +| 

Pj_ 

M 

P*> 

U-xJ 

Fluids:  R12,  R22.  R113 
7000  <  Refe  <  53000 


8.42 


Shah  (1979) 


KD» 


Nu} 


0.023  Re*  Pr°‘ 


,,  ,0.0 

('-■*)  + 


3.8.x0'76  (l-x)° 


D  =  7-40  mm 

Fluids:  water.  R1 1.  R12, 

R22, 

R 1 1 3,  methanol,  ethanol, 
benzene,  toluene, 
trichloroethylene _ 


13.32 


Moser  el  al. 
(1998) 


KD!'  _  Nu, 
kf  Nu,  j 


0.0994“ 1  Re/  Re],*08756'  Pr" 815 


(l.581nRet.9-3.28)(2.581nRea/+13.7Pr;  3-19.l) 
C,  =0.126Pr/448  ,  C2  =-0.1 13Pr;“5W,  Kteq=^Friedd Re* 


£>  =  3.14-20  mm 
Fluids:  R1 1,  R12,  R125, 
R22, 

R 1 34a,  R4 1 0A 
1 197  data  points _ 
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Dobson  and 
Chato 
(1998) 
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X..  = 


Nu, 
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2.22 

^  V-O.SV 
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> 

,0.9 

f  ^0 
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—  l 

U, 

l  X  J  1 
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Wang  et  al. 
(2002) 


KD> •  _  f  Nu, 
Nu, 


0.0274  Pr,  Re,  x 


0.6792  0.2200  ^R 


X.. 


£>/,  =  1 .46  mm 
Fluid:  R134a 
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flS2  =1.376 +  8X1;665 

rR  II 
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Koyama  et 
al. 

(2003) 


Mi  =  ^j0.0152(l  +  0.6Pr;8)-^Re7 
/g  =  l  +  2l[l-  exp(-0.3 1 90, ,  )]*„  +  X2„ 
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Huang  et  al. 

(2010) 
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Figure  3.14  Comparison  of  present  experimentally-determined  average  FC-72 
condensation  heat  transfer  coefficient  data  with  predictions  of  annular  flow 
macro-channel  condensation  correlations. 
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Figure  3.15  Comparison  of  present  experimentally-determined  average  FC-72 
condensation  heat  transfer  coefficient  data  with  predictions  of  annular  flow 
mini/micro-channel  condensation  correlations. 
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Notably,  the  correlations  of  Cavallini  and  Zecchin  (1974),  Shah  (1979),  and  Dobson  and 
Chato  (1998)  predict  the  present  FC-72  data  with  MAEs  of  8.42%,  13.32%,  and  10.65%, 
respectively. 

Figure  3.15  compares  the  present  FC-72  average  condensation  heat  transfer 
coefficient  data  with  predictions  of  annular  condensation  heat  transfer  correlations 
intended  for  mini/micro-channels.  Wang  et  al.  (2002)  developed  a  boundary  layer-based 
correlation  for  annular  flow  in  which  the  two-phase  frictional  multiplier  and 
dimensionless  boundary  layer  temperature  were  evaluated  for  R134a  condensing  inside  a 
horizontal  rectangular  multi-channel  with  a  hydraulic  diameter  of  1.46  mm.  The 
correlation  of  Koyama  et  al.  (2003)  is  based  on  Haraguchi  et  al.  (1994b)’s  correlation  for 
annular  flow,  which  is  derived  from  boundary  layer  analysis  for  condensation  of  R22, 
R134a,  and  R123  inside  a  horizontal  8.4  mm  diameter  circular  tube  and  Mishima  and 
Hibiki’s  (1996)  correlation  for  the  two-phase  frictional  multiplier.  Huang  et  al.'s  (2010) 
correlation  is  based  on  Haraguchi  et  al.' s  (1994a)  frictional  multiplier  correlation  for  a 
diameter  of  8.4  mm  and  Haraguchi  et  al.' s  (1994b)  annular  condensation  correlation, 
which  Huang  et  al.  modified  with  their  own  experimental  data.  Although  the  correlations 
in  Fig.  3.15  are  intended  for  mini/micro-channels,  they  show  relatively  poor  agreement 
with  the  present  FC-72  data  compared  to  the  top  four  macro-channel  correlations  in  Fig. 
3.14. 


3.3.4.  New  Heat  Transfer  Correlation 


For  a  shear-dominated  annular  condensing  flow,  the  local  condensation  heat 
transfer  coefficient  of  liquid  film  can  be  obtained  from  the  following  relation, 


h 


<7„ 


lp,cir 


T  -T 

sal  w 


P/cpj1 

t; 


(3.28) 


where  the  friction  velocity  and  dimensionless  boundary  layer  temperature  are  defined, 
respectively,  as 

«*  =  JT*/Pf  >  (3.29a) 


and 


77  =  f’  S- 
S  Jo  ql 


r  ,  .  \ 

1  1  £ 
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v  Pry  Prr  v 


f  ) 


dy+  . 


(3.29b) 
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The  wall  shear  stress  can  be  expressed  as  a  function  of  the  single-phase  frictional 
pressure  drop  based  on  liquid  flow  and  the  two-phase  frictional  multiplier, 


'dP' 


\dz  j 


D, 


dP' 

dz 


D 


2 

f  ■ 


(3.30) 


The  single-phase  frictional  pressure  drop  based  on  liquid  flow  can  be  expressed  as 

2/G2(  \-x)2of 


dP 

dz 


D„ 


(3.31) 


The  correlation  of  Mishima  and  Hibiki  (1996),  which  is  based  on  adiabatic  air- 
water  two-phase  flow  inside  1-4  mm  diameter  circular  tubes,  showed  good  predictions  of 
experimental  pressure  drop  data  in  small  channels  such  as  circular  tubes  (Cavallini  et  al., 
2005),  trapezoidal  channels  (Quan  et  al.,  2008),  and  rectangular  channels  (Koyama  et  al., 
2003;  Qu  and  Mudawar,  2003a;  Park  and  Hrnjak,  2009).  Mishima  and  Hibiki’s 
correlation  also  provided  the  best  predictions  of  the  present  FC-72  pressure  drop  data 
(MAE  of  27.15%)  among  all  the  separated  flow  correlations  examined  in  Chapter  3.2. 
Therefore,  the  two-phase  frictional  multiplier  in  Eq.  (3.30)  is  based  on  Mishima  and 
Hibiki’s  correlation 


,  2l[l-exp(-0.319Z>,)] 

Tf  V 


Dh[mm \ . 


(3.32) 


It  was  shown  earlier  that  most  macro-channel  correlations  intended  for  circular 
tubes  provide  good  agreement  with  the  present  square  micro-channel  condensation  heat 
transfer  data.  This  suggests  that,  unlike  non-annular  flow  regimes,  there  are  no  major 
differences  in  condensation  behavior  between  small  and  large  channels  in  the  annular 
regime.  Furthermore,  the  low  surface  tension  of  FC-72  tends  to  maintain  nearly  uniform 
liquid  film  thickness  along  the  inner  walls  of  the  square  channel,  barely  affected  by 
corner  effects.  Therefore,  the  following  relation  by  Haraguchi  et  al.  (1994b)  for 
condensation  of  R22,  R134a,  and  R123  inside  an  8.4-mm  horizontal  circular  tube  with  G 
=  200  -  400  kg/m2s  is  used  for  determination  of  the  dimensionless  boundary  layer 
temperature  in  the  present  square  micro-channel. 


Ts+  =  (0. 1  +  0.06Pr“'8 )"'  Pr,  Re0/3 .  (3.33) 

Substituting  Eqs.  (3.29)  -  (3.33)  into  Eq.  (3.28),  and  accounting  for  the  effects  of 
three-sided  cooling  via  Eq.  (3.26),  the  local  heat  transfer  coefficient  for  annular 
condensation  with  three-sided  cooling  can  be  expressed  as 
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^=fe)0-1  +  0-MPr“)Pr;lRc7l,3c^ 

(3.34) 

where  X  is  the  Martinelli  parameter.  The  Fanning  friction  factor  based  on  liquid  flow  can 
be  determined  from  (Shah  and  London,  1978;  Incropera  and  Dewitt,  2002). 

/Re/  =  24(1-1.3553/?+ 1.9467/?2  -1.70 12 /?3  +  0.9564/?4  -0.2537/?5) 

for  Re  f  <  2,000 ,  (3.35a) 

f  =  0.079 Re/25  for  2,000  <  Re,  <  20,000 ,  (3.35b) 

/  =  0.046 Re/2  for  Re,  >  20,000,  and  Re,  =  G(\-x)Dh//uf  .  (3.35c) 

Figure  3.16  shows  the  new  condensation  heat  transfer  correlation  accurately 
predicts  the  present  FC-72  average  heat  transfer  data  with  a  mean  absolute  error  of  9.0%, 
although  some  low  mass  velocity  data  include  non-annular  data. 

Figure  3.17  compares  the  present  local  heat  transfer  coefficient  data  for  G  =  306 
kg/m2s  and  mw  =  5  g/s  with  predictions  of  the  new  correlation,  Eq.  (3.34),  and  previous 
correlations  given  in  Table  3.4.  The  new  correlation  shows  best  predictions  of  the 
present  data  both  in  magnitude  and  trend.  The  correlations  of  Akers  et  al.  (1958),  Wang 
et  al.  (2002)  and  Huang  et  al.  (2010)  over-predict  the  data  by  an  appreciable  margin. 
While  the  Shah  (1979)  correlation  predicts  an  unrealistic  trend  of  decreasing  heat  transfer 
coefficient  in  the  high  vapor  quality  region  where  the  film  is  thinnest,  this  correlation 
provides  good  overall  prediction  of  the  average  heat  transfer  coefficient,  evidenced  by  a 
MAE  of  13.32%. 


In  order  to  further  assess  the  accuracy  of  the  new  annular  condensation  heat 
transfer  correlation,  a  total  923  data  points  for  condensation  in  mini/micro-channels  were 
amassed  from  eight  sources  which  are  described  in  Table  3.5.  Figure  3.18  shows  all  923 
data  points  in  a  plot  of  Martinelli  parameter  based  on  turbulent  liquid-turbulent  vapor 
versus  modified  Weber  number,  which  is  defined  by  Soliman  (1986)  as 


We  =  2.45 


S«f(l  +  L09XrT 


for  Re,  <  1250, 


(3.36a) 
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where  the  Suratman  number,  Sug,  and  the  turbulent-turbulent  Martinelli  parameter,  X„, 
are  defined  as 


Sug 


PS°D 


(3.37a) 


and 


*„  = 


(  Y 

°-1  f  ,0.9 

f  \ 

V 

l~X) 

0.5 


(3.37b) 


respectively  (Soliman,  1986).  The  vapor  and  liquid  Reynolds  numbers  of  the  923  data 
points  are  also  symbolized  in  Fig.  3.18  with  a  transition  Reynolds  number  value  of  2000. 
Although  all  the  data  are  plotted  versus  the  turbulent-turbulent  Martinelli  parameter,  they 
are  classified  clearly  into  the  following  four  different  zones  by  two  fitted  lines  of 
We*  =3  and  We *  =6.3X/068  , 


Ref  <  2000  and  Reg  <  2000:  We * 

<  3  and  We*  <  6.3X,,'0'68 . 

(3.38a) 

Ref  <  2000  and  Rea  >  2000: 

3  <We  <  6.3X,,'0'68 . 

(3.38b) 

Re,  >  2000  and  Rea  <  2000: 

6. 3 A',,’0'68  <We*  <3. 

(3.38c) 

Re,  >  2000  and  Rea  >  2000:  We 

>3  and  We  >6.3A'/068  . 

(3.38d) 

Table  3.5  Database  for  two-phase  condensation  heat  transfer  coefficient  in 

mini/micro-channels. 


Author(s) 

Dh  [mm] 

Fluid(s) 

Geometry* 

All  data  / 

annular  data 

Hirofumi  and  Webb  (1995) 

0.96-2.13 

R134a 

C/R  multi 

62/61 

Zhang  (1998) 

2.13,3.25 

R134a,  R22,  R404A 

C  single/multi 

66  /  65 

Yan  and  Lin  (1999) 

2.0 

R134a 

C  multi 

78/60 

Wang  et  al.  (2002) 

1.46 

R134a 

R  multi 

410/268 

Cavallini  et  al.  (2005) 

1.4 

R134a,  R410A 

R  multi 

59/56 

Park  and  Hmjak  (2009) 

0.89 

C02 

C  multi 

52/37 

Matkovic  et  al.  (2009) 

0.96 

R 1 34a,  R32 

C  single 

161  / 131 

Huang  et  al.  (2010) 

1.6,  4.18 

R410A 

C  single 

35/35 

*  C:  circular,  R:  rectangular,  'annular  data  (smooth-annular,  wavy-annular  and  transition)  correspond  to  We’  >  7X„02. 
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Figure  3.16  Comparison  of  present  experimentally-determined  average  FC-72 
condensation  heat  transfer  coefficient  data  with  predictions  of 
the  new  annular  condensation  correlation. 
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Figure  3.17  Variation  of  present  experimentally-determined  local  FC-72  condensation 
heat  transfer  coefficient  data  with  thermodynamic  equilibrium  quality  for 
G  =  306  kg/  m2s  and  mw  =  5  g/s  compared  to  predictions  of  new  correlation 
and  previous  annular  condensation  correlations. 
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Figure  3.19  Proposed  flow  regime  map  based  on  present  FC-72  flow  visualization 
experiments  in  square  micro-channels  with  Dh  =  1  mm  plotted  alongside  639 
mini/micro-channel  data  points  corresponding  to  Jg*  >  2.5  from  eight  sources. 
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Based  on  the  FC-72  flow  visualization  data  presented  in  Chapter  3.1,  the 
following  lines  are  fitted  for  boundaries  between  different  flow  regimes: 

Smooth-annular  to  wavy-annular:  We’=90X„05.  (3.39a) 

Wavy-annular  to  transition:  We*  =  24X,0Al.  (3.39b) 

Transition  to  slug:  We *  =  IX,,02 .  (3.39c) 

Cavallini  et  al.  (2002)  proposed  that  the  flow  is  annular  for  Jg  >  2.5,  and 
intermittent  and  slug  for  J*  <  2.5,  based  on  previous  flow  regime  maps,  where  the 
dimensionless  superficial  vapor  velocity  is  defined  as 

j',=g*NpM-pXd-  <3-4°) 

Figure  3.19  shows  639  mini/micro-channel  annular  data  points  from  the  eight 
aforementioned  sources  excluding  all  data  with  dimensionless  superficial  gas  velocities 
below  2.5.  Although  a  few  data  points  do  appear  for  We*  <  7X„o:,  the  boundary  line 
between  the  transition  and  slug  regimes  from  the  present  study  is  in  general  agreement 
with  that  from  Cavallini  et  al.. 


Table  3.6  Comparison  of  mini/micro-channel  database  for  annular  flow  (smooth-annular, 
wavy-annular  and  transition)  with  predictions  of  annular  condensation 
heat  transfer  correlations. 


Authors) 

Mean  absolute  error  [%] 

Akers 
et  al. 
(1958) 

Cavallini 

and  Zecchin 
(1974) 

Shah 

(1979) 

Moser 
et  al. 
(1998) 

Dobson 

and  Chato 
(1998) 

Wang 
et  al. 
(2002) 

Koyama 
et  al. 
(2003) 

Huang 
et  al. 
(2010) 

New 

correl. 

Hirofutni  and  Webb  ( 1 995) 

89.1 

84.9 

72.4 

39.1 

77.9 

12.3 

30.2 

30.0 

29.3 

Zhang  (1998) 

27.6 

32.6 

22.5 

10.6 

27.8 

33.7 

50.4 

36.6 

15.5 

Yanand  Lin  (1999) 

84.5 

16.8 

12.0 

11.8 

16.4 

21.0 

48.3 

39.6 

13.4 

Wang  et  al.  (2002) 

27.2 

59.8 

45.8 

22.2 

46.5 

14.1 

45.5 

38.4 

16.8 

Cavallini  et  al.  (2005) 

53.3 

33.5 

23.7 

12.6 

28.7 

33.6 

49.9 

28.8 

17.1 

Park  and  Hmjak  (2009) 

178.3 

73.6 

54.6 

40.8 

67.7 

14.7 

30.0 

19.7 

11.0 

Matkovic  et  al.  (2009) 

59.6 

29.9 

18.8 

14.0 

26.8 

34.5 

48.4 

22.4 

21.6 

Huang  et  al.  (2010) 

12.8 

19.1 

15.2 

15.8 

17.4 

48.2 

57.3 

58.9 

20.8 

Total 

90.1 

47.0 

35.4 

20.1 

39.6 

23.3 

45.5 

34.0 

18.3 
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Figure  3.20  Comparison  of  predictions  of  new  annular  condensation  heat  transfer 
correlation  with  prior  experimental  mini/micro-channel  data. 
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In  order  to  isolate  data  corresponding  to  flow  regimes  with  a  clearly  identifiable 
annular  liquid  film  (smooth-annular,  wavy-annular  and  transition)  to  assess  the  accuracy 
of  the  new  correlation,  all  210  data  points  (from  a  total  of  923  data  pints)  that  correspond 
to  We*  <  7X„tu  are  excluded.  Table  3.6  shows  both  the  total  number  of  data  points  from 
each  of  the  eight  previous  sources  as  well  as  corresponding  number  of  data  points  that  are 
deemed  to  fall  into  the  present  broad  classification  of  annular  flows. 

Table  3.6  compares  the  predictive  accuracy  of  previous  correlations  and  the  new 
annular  flow  correlation  against  the  mini/micro-channel  data  that  fall  into  the  broad 
classification  of  annular  flows.  Although  the  correlations  of  Moser  et  al.  (1998)  and 
Wang  et  al.  (2002)  provide  good  predictive  capability  against  the  entire  annular  database, 
with  MAEs  of  20.1%  and  23.3%,  respectively,  they  are  far  less  accurate  in  predicting  the 
present  FC-72  data,  with  corresponding  MAEs  of  31.8%  and  88.0%,  respectively.  The 
new  annular  heat  transfer  correlation  accurately  predicts  each  of  the  eight  previous 
sources  with  a  MAE  less  than  30%. 

Figure  3.20  shows  the  new  annular  correlation  provides  the  best  prediction  of  713 
annular  mini/micro-channels  data  points  with  an  overall  MAE  of  18.3%  (17.6%  for 
multi-channels  and  19.8%  for  single  channels).  Most  data  points  are  shown  falling 
within  ±30%  of  predicted  values. 
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CHAPTER  4.  THEORETICAL  MODEL 


4.1.  Model  Development 


4.1.1.  Key  Findings  from  Flow  Visualization  Experiments 

Qu  and  Mudawar  (2003a)  showed  that  droplet  entrainment  and  deposition  play  a 
crucial  role  in  mass,  momentum  and  heat  transfer  in  annular  evaporating  micro-channel 
flows.  As  shown  in  Fig.  3.8(a),  droplets  are  formed  by  shattering  of  liquid  from  the 
micro-channel’s  upstream.  However,  the  flow  visualization  results  from  the  present 
authors’  study  have  shown  that  no  droplets  are  formed  in  the  annular  region  of  micro- 
channel  condensing  flow  (refer  to  Chapter  3.1).  Therefore,  these  effects  are  neglected  in 
the  development  of  the  present  theoretical  annular  flow  model. 

The  authors’  recent  micro-channel  condensation  studies  also  showed  that  the 
annular  region  consists  of  two  distinct  sub-regimes.  The  upstream,  smooth-annular 
regime  features  a  thin  smooth  liquid  film.  However,  the  annular  film  in  the  wavy-annular 
regime  is  marred  by  interfacial  waves.  Nonetheless,  these  waves  were  fairly  sinusoidal  in 
profile  and  do  not  shatter  or  produce  entrained  droplets,  nor  do  they  evolve  into  large 
waves  as  is  commonly  observed  in  annular  macro-channel  flows. 


4.1.2.  Model  Assumptions 

A  key  challenge  in  developing  the  theoretical  model  is  determining  the 
distribution  of  liquid  along  the  channel’s  circumference.  Gravity  effects  are  negligible  in 
micro-channels  compared  to  macro-channels,  given  the  large  shear  stresses  encountered 
in  the  former.  This  precludes  any  preferential  accumulation  of  liquid  based  on 
orientation  of  the  flow  relative  to  gravity.  However,  the  surface  tension  of  the  working 
fluid  can  have  a  significant  influence  on  the  liquid  distribution.  Figure  3.8(b)  contrasts 
the  liquid  distribution  in  a  rectangular  micro-channel  for  the  case  of  a  very  low 
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surface  tension  fluid  such  as  FC-72  (0.0084  N/m  at  1  bar)  compared  to  that  of  a  relatively 
high  surface  tension  fluid  such  as  water  (0.059  N/m  at  1  bar).  Notice  how  low  surface 
tension  fluids,  which  include  most  dielectric  fluids  used  in  electronic  cooling  applications, 
tend  to  maintain  a  nearly  uniform  liquid  film  thickness  around  the  channel  circumference. 
The  present  study  will  be  based  on  this  assumption  since  it  concerns  mostly  low  surface 
tension  electronic  cooling  fluids. 

Figure  4.1  provides  a  schematic  representation  of  liquid  film  condensation  in  a 
rectangular  micro-channel  with  three-sided  wall  cooling  and  relevant  nomenclature  used 
in  the  development  of  the  annular  condensation  model.  Following  are  key  assumptions 
of  the  model: 

(1)  The  annular  condensing  flow  is  steady,  incompressible,  and  concurrent. 

(2)  Gravitational  effects  are  negligible. 

(3)  Both  entrainment  and  deposition  of  liquid  droplets  are  negligible. 

(4)  Pressure  is  uniform  across  the  micro-channel’s  cross-sectional  area. 

(5)  Thermophysical  properties  are  based  on  local  saturation  pressure. 

(6)  Mean  vapor  velocity  is  assumed  across  the  vapor  core’s  cross-sectional  area. 

(7)  Mass  transfer  occurs  only  at  the  interface  between  the  vapor  core  and  the  liquid  film. 

(8)  The  liquid  film  interface  is  smooth. 

(9)  Liquid  film  thickness  is  uniform  around  the  channel’s  circumference. 

(10)  Axial  momentum  changes  in  the  liquid  film  are  negligible. 

(11)  The  liquid  film  flow  can  be  laminar  or  turbulent. 

(12)  The  circumferential  heat  flux  is  uniform  along  the  vapor-liquid  interface. 

The  boundary  condition  of  base  heat  flux,  q"base,  beneath  the  micro-channels  is 
determined  from  local  temperature  measurements  obtained  by  arrays  of  thermocouples 
embedded  in  two  separate  planes  in  the  copper  block.  The  local  temperatures  in  the 
copper  block  are  fitted  to  a  third-  or  fourth-order  polynomial  function  of  stream-wise 
distance,  and  the  corresponding  variations  of  micro-channel  base  heat  flux,  q"ba*e,  and 
base  temperature,  Thaxe,  are  obtained  by  assuming  one-dimensional  heat  conduction  along 
the  direction  perpendicular  to  the  top  surface  of  the  copper  block  (Kim  and  Mudawar, 
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Condensing  film 


Figure  4. 1  Schematic  representations  of  (a)  condenser  with  rectangular  micro-channels 
and  (b)  liquid  film  condensation  in  micro-channel  with  three-sided  wall  cooling. 
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2010a).  In  the  superheated  vapor  region  near  the  channel  inlet,  sensible  heat  loss  of  FC- 
72  can  be  calculated  from  the  following  energy  balance, 


?L,fc  +  »',K  =  ^c,,,Ar/.  (4.1) 

The  onset  of  two-phase  condensation  is  calculated  to  occur  at  the  distance  from 
the  channel  inlet  where  the  fluid  temperature  reaches  the  saturation  temperature 
corresponding  to  the  local  pressure.  The  quality  in  the  superheated  region  can  be 
determined  from 


jc  =  1 


cp.g  (^/  Tsll, ) 


h 


(4.2) 


fg 


where  7}  and  Tsa,  are  the  local  bulk  fluid  and  saturation  temperatures,  respectively. 


4.1.3.  Control  Volume  Analysis 


To  develop  a  theoretical  control-volume-based  model,  mass  and  momentum 
conservation  are  first  applied  to  control  volumes  encompassing  the  liquid  film  and  the 
vapor  core  separately.  For  the  two-phase  condensing  region,  mass  conservation  for  the 
liquid  film  and  the  vapor  core  can  be  expressed,  respectively,  as 


d  m , 
dz 


=  0, 


(4.3) 


and 


d  nr , 
dz 


+  r/g-°' 


(4.4) 


where  the  mass  flow  rates  of  the  liquid  film  and  the  vapor  core,  and  the  rate  of  mass 
transfer  due  to  condensation  are  defined,  respectively,  as 


and 


>iif  =  2 pf  £  uf  [(Hch  -  2 y)  +  (Wch  -  2y)\iy , 
mg=psus(Wch-2S\Hch-28), 

*:(2  Hch+Wth)  qlse(Wch+Ws ) 


(4.5) 

(4.6) 

(4.7) 

‘fg  "fg 

Applying  momentum  conservation  to  the  liquid  film  element  illustrated  in  Fig. 
4.2(a)  yields 


r/g  = 


h, 


h . 


r 


~r fHuiAz  =  PAf.- 


dP 


P  +  —  Az 
V  dz 


A  f,»  —  TPf,y  Az  +  T,  Pf'g&Z, 


(4.8) 
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where  the  flow  area.  A/;*,  local  perimeter,  Pfy,  and  interfacial  perimeter,  P/;s,  can  be 
expressed,  respectively,  as 

Af,  =  -lytw*  -2  y)-(H«-2S\W« -25),  (4.9) 

Pf,y=  2[(W«, -2v)+(Wch  -2y)],  (4.10) 

and  P/?e=2[(Hct-2S)+{Wci-25)].  (4.11) 

Rearranging  Eq.  (4.8)  gives 


r  ,r,^Ar.  T.Prx+r,U 


T= 


dP_ 

dz 


V  J 


p, 


/,*  +_i*  f.8  '  '/s“i 


Pr 


(4.12) 


f.y  1  f.y 

Allowing  for  turbulence  in  the  condensing  film,  the  local  shear  stress  in  the  film 
can  be  expressed  as 

r  \ 


T=M, 


1 


v  v/. 


du . 


dy 


(4.13) 


where  em  is  the  eddy  momentum  diffusivity.  Substituting  Eq.  (4.13)  into  Eq.  (4.12)  and 
integrating  yield  the  velocity  profile  across  the  liquid  film. 


uf{y)  = 


-f 

dP'' 

if y/s  Af,> 

(  ^ 

1  +  — 

-1 
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ry] 

Mf  V 

dz  y 

Jo  P 

‘f.y 

l  VfJ 

kS) 

Pf 


i 

-1 

d\ 

ry) 

pf. 

K  VfJ 

\S) 

(4.14) 


The  interfacial  velocity  can  be  determined  by  setting  y  =  <$>’ in  Eq.  (4. 14). 


u,  -- 


f  dP)fAfs 

f  > 

1  +  — 

S'  "N 
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+  Tipf.st 

r 

f  1 
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v<?y 


(4.15) 


Substituting  Eq.  (4.14)  into  Eq.  (4.5)  yields  the  following  relation  for  pressure 
gradient. 


P,mf 

dP _  p<  5~ 

dz 


(r/  Pf,5  +  r/g Ui  )J0' 
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■y/S  1 
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'  _  V1 
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v  V(  J 
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(4.16) 


The  double  integral  terms  in  the  above  equation  can  be  expanded  as  follows, 
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(4.17) 


Applying  momentum  conservation  to  a  rectangular  element  of  the  vapor  core 
illustrated  in  Fig.  4.2(b)  yields 


P*ulA* 


P\~ 


+  diyPs2A^  PgU*Ag  +  T*U,6Z  = 

PAl+^ Az 
*  dz 


(  1  dP 

—T:  Pr  ,AZ  + 

P  + - Az 

t  / 

V  2  dz 

dA 


(4.18) 


dz 


LAz 


Rearranging  Eq.  (4.18)  and  neglecting  the  second  order  term  of  A z  yield  the  following 
relation  for  interfacial  shear  stress  between  the  condensing  film  and  vapor  core, 

dP'\  4pXAJ 


1 


r.  =  ■ 


f.s 


A. 


dz 


s  8  *’—T  u 

1  fSUi 


dz 


(4.19) 


where  the  flow  area  of  the  vapor  core  is  A„  =  (Wch  -  28 \Hch  -28'). 

The  interfacial  shear  stress  is  the  result  of  velocity  differences  between  the  vapor 
core  and  interface,  modified  by  the  influence  of  interfacial  momentum  transfer  due  to 
condensation  along  the  vapor-liquid  interface  due  to  condensation;  the  later  is  obtained 
using  a  treatment  by  Wallis  (1969), 


Ti  =  f 
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Pg  \~‘g 


k  ~u4 
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(w 


fg 


2  P. 


(4.20) 


f.s 


The  interfacial  friction  factor,  f,  can  be  determined  from  relations  by  Shah  and  London 
(1978), 

f\ Ret.  =24(l-1.3553^.  +  1.9467/?/  -1.7012#.3  +0.9564#. 4  -0.2537/?/) 

for  Re(.  <2,000,  (4.21a) 

f  =  0.079Re/25  for  2,000  <  Re(.  <  20,000 ,  (4.2 1  b) 

and  f  =  0.046 Re/2  for  Re(.  >  20,000 ,  (4.21c) 

where  /?c,  Rec,  and  Dhc  are  the  aspect  ratio,  effective  Reynolds  number,  and  hydraulic 
diameter  of  the  vapor  core,  respectively,  which  are  given  by 

fV.,-28 
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and 


Re 


H 


(4.23) 

(4.24) 


4.1.4.  Turbulence  Model 


Equation  (4.13)  can  be  expressed  in  nondimensional  form  as 
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(4.25) 

(4.26) 

(4.27) 

(4-28) 


Based  on  the  Prandtl  mixing  length  theory,  the  eddy  diffusivity  can  be  expressed 
in  terms  of  the  turbulent  mixing  length  according  to  the  relation 

.2  du  + 


=  r 


(4.29) 


vf  dy 

A  turbulent  mixing  length  relation  originally  proposed  by  Van  Driest  (1956)  was 
modified  by  Kays  (1972,  1980)  to  the  following  form. 


P  -  Ky 


1  -  exp 


_LZ_ 

1  T...  A H 


(4.30) 


where  the  Von-Karman  constant  of  K  =  0.4  is  used  for  all  subsequent  calculations,  and 
the  constant  A*  is  given  by  Kays  (1972,  1980), 


A +  =  26 


-0.5  _  -1.5 


l  +  30AS/jfp~  tw 


dP 

dz 


(4.31) 


Substituting  Eqs.  (4.25)  and  (4.30)  into  Eq.  (4.29)  yields  the  following  eddy  diffusivity 
profile. 
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1  +  4  K2y+2 


1  -  exp 


(4.32) 


In  order  to  account  for  interfacial  dampening  in  the  condensation  film  due  to 
surface  tension  suppression  of  turbulent  eddies,  a  dampening  term,  (l -y+/S+)" ,  is 
included  in  the  above  eddy  diffusivity  profile,  where  the  influence  of  the  parameter  n  will 
be  discussed  later.  The  complete  form  of  eddy  momentum  diffusivity  distribution  in  the 
shear-driven  film  is  expressed  as 
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where,  based  on  Eq.  (4.12), 
r  ^2(Hcl+Wj 
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(4.34) 


4.1.5.  Determination  of  Heat  Transfer  Coefficient 


Heat  flux  across  the  liquid  film  is  related  to  the  liquid  temperature  gradient  by  the 
relation, 
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where  T *  is  the  dimensionless  temperature  defined  as 


(4.35) 


T+ 


PfcP,fu‘(T-T») 


(4.36) 


and  Prr  is  the  turbulent  Prandtl  number  ( eJXh ),  which,  as  discussed  by  Mudawar  and  El- 
Masri  (1986),  can  be  evaluated  from  the  experimental  data  of  Ueda  et  al.  (1977), 


Prr  =  1.4exp 


5 + 


/ 


0.66. 


(4.37) 


Integrating  Eq.  (4.35),  the  dimensionless  temperature  profile  across  the  liquid  film  can  be 
expressed  as 
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r=f 
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1  em 
Prr  VfJ 


dy 


(4.38) 


Under  three-sided  cooling  conditions  for  trapezoidal  micro-channels,  Quan  et  al. 
(2010)  developed  a  semi-analytical  model  for  the  annular  condensation  heat  transfer 
coefficient  based  on  the  assumption  of  uniform  circumferential  heat  flux  in  the  vapor- 
liquid  interface.  Adopting  a  heat  flux  relation  similar  to  that  of  Quan  et  al.,  the  present 
model  for  condensation  in  rectangular  micro-channels  with  three-sided  cooling  walls 
employs  the  following  energy  balance. 
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Therefore,  the  local  condensation  heat  transfer  coefficient  can  be  expressed  as 


(4.39) 
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4.1.6.  Calculation  Procedure 

The  model  equations  presented  in  the  previous  sections  are  solved  numerically 
using  a  finite  difference  technique.  The  axial  distance  is  divided  into  small  zlz  increments 
and  calculations  are  repeated  staring  at  the  upstream  location  where  the  film  is  initiated. 
The  calculation  procedure  is  as  follows: 

1.  smj vf  and  5  are  both  set  to  zero  at  the  film’s  upstream  location. 

2.  An  initial  guess  of  emj  vf  at  the  next  Az  location  is  made  with  the  values  of  the  wall 
shear  stress,  rw,  interfacial  shear,  r,-,  interfacial  velocity,  «„  and  pressure  gradient,  -dPIdz, 
at  the  node  immediately  upstream  using  Eq.  (4.33). 

3.  An  initial  value  of  the  film  thickness  8  is  assumed. 

4.  The  interfacial  velocity,  m„  interfacial  friction  factor,  f,  interfacial  shear,  r,-,  and 
pressure  gradient,  - dPIdz ,  are  calculated  using  Eqs.  (4.15),  (4.21),  (4.20),  and  (4.16), 
respectively. 

5.  Convergence  is  checked  by  comparing  the  two  sides  of  Eq.  (4.19).  If  the  sides  are  not 
equal,  sm/vf  is  updated  with  new  values  of  tw,  t,,  «„  and  -dPIdz  at  the  present  node,  and 
steps  (3)-(5)  are  repeated  until  the  correct  value  for  t>is  found. 
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6.  The  values  of  emj vf ,  «+,  A,  and  h,p  are  calculated  using  Eqs.  (4.33),  (4.14),  (4.38), 
and  (4.40),  respectively.  Steps  (2)-(6)  are  then  repeated  for  the  next  downstream  node. 
The  procedure  is  continued  until  the  last  node  is  reached. 


4.1.7.  Simplified  Model 


Since  the  thickness  of  the  condensing  film  is  generally  small  compared  to  the 
micro-channel  hydraulic  diameter,  a  simplified  model  can  be  derived  using  the  following 
assumption, 

P,.y=P,J=Pf='AHct+Wct).  (4.41) 

Then,  the  flow  area  and  mass  flow  rate  of  the  liquid  film  can  be  expressed,  respectively, 
as 


Af,=pAs-y)> 

mf  =PfPf[)  ufdy . 


(4.42) 


and  m ,  -  p , P,  \  u,dv.  (4.43) 

Applying  momentum  conservation  to  a  rectangular  annular  element  of  the  liquid 
film  yields  the  following  shear  stress  relation, 
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Similar  to  the  procedure  used  in  previous  section,  the  velocity  profile,  interfacial 
velocity,  and  pressure  gradient  can  be  simplified,  respectively,  as 
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As  discussed  later,  the  predictive  differences  of  pressure  drop  average  only  1.8%, 
with  a  maximum  difference  of  3.8%;  corresponding  differences  for  the  average  heat 
transfer  coefficient  are  0.5%  and  1.0%,  respectively.  Although  both  models  yield  very 
close  results,  all  of  the  subsequent  calculations  are  based  on  the  complete  model  unless 
indicated  otherwise. 


4.2.  Model  Results 


4.2.1.  Effect  of  Turbulent  Dampening  Term 


Mudawar  and  El-Masri  (1986)  examined  the  interfacial  dampening  effects  on  the 
turbulence  mixing  length  for  free-falling  liquid  films.  They  derived  the  following 
turbulence  mixing  length  profile  based  on  an  eddy-diffusivity  profile  measured 
experimentally  by  Ueda  et  al.  (1977)  for  open-channel  flows, 
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Substituting  Eqs.  (4.25)  and  (4.48)  into  Eq.  (4.29)  gives  the  following  eddy  diffusivity 
profile. 
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It  should  be  noted  that,  since  the  mixing  length  profile  of  Mudawar  and  El-Masri  is  valid 
only  for  smjvf>  1  and  negligible  interfacial  shear  stress,  their  eddy  diffusivity  profile 
cannot  be  used  for  micro-channel  condensation,  where  turbulence  effects  are  relatively 
weak  and  interfacial  shear  is  the  primary  driving  force  for  film  motion. 

The  present  eddy  diffusivity  profile  given  by  Eq.  (4.33)  yielded  values  for 
dimensionless  film  thickness,  y+  =  8r,  up  to  about  16  at  the  micro-channel  exit,  which  is 
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close  to  the  outer  edge  of  the  viscous  sublayer  for  typical  turbulent  flows.  The  present 
diffusivity  profile  is  justified  by  the  fact  that  Van-Driest  models  are  applicable  near  the 
solid  wall  region  where  y+  <  30.  The  dampening  coefficient,  (l  -  y+  /S* )" ,  in  Eq.  (4.33) 
accounts  for  the  vanishing  turbulence  effects,  if  any,  at  the  film  interface.  Without  this 
term,  the  conventional  Van-Driest  profile  cannot  account  for  the  interfacial  dampening. 

Using  Eq.  (4.33),  four  different  powers  of  the  dampening  term,  (l  -  y*  /<?+)" ,  n  = 
0  (corresponding  to  zero  dampening  as  expressed  by  Eq.  (4.32)),  0.1,  0.5,  and  1.0,  are 
examined.  Figures  4.3(a),  4.3(b)  and  4.3(c)  show  the  impact  of  the  dampening  term  on 
the  distributions  of  eddy  diffusivity,  velocity,  and  temperature,  respectively,  across  the 
FC-72  film  for  G  =  367  kg/m2s  and  Wlw  =  6  g/s.  Figure  4.3(a)  shows  that,  without  the 
dampening  term,  the  eddy  diffusivity  increases  monotonically  with  increasing  y+,  while 
with  the  dampening  term,  the  eddy  diffusivity  is  reduced  to  zero  at  the  film  interface. 
Increasing  the  magnitude  of  the  exponent  n  causes  the  effect  of  interfacial  dampening  to 
penetrate  deeper  towards  the  wall.  Figure  4.3(b)  shows  that  increasing  n  increases  the 
liquid  velocity  near  the  interface.  Figure  4.3(c)  shows  a  similar  trend  of  increasing  liquid 
temperature  near  the  interface  with  increasing  n.  An  interesting  aspect  of  the  temperature 
profile  is  the  unusual  temperature  rise  in  the  immediate  vicinity  of  the  film  interface 
caused  by  the  eddy  diffusivity  term  in  Eq.  (4.35)  approaching  zero. 

Although  Fig.  4.4(a)  shows  that  different  n  values  yield  fairly  similar  film 
thickness  results,  Fig.  4.4(b)  shows  that  increasing  the  power  increases  the  temperature 
gradient  near  the  interface,  which  can  produce  a  measureable  decrease  in  the  heat  transfer 
coefficient.  An  optimum  value  of  n  =  0. 1  was  ascertained  by  comparing  predictions  of 
the  local  as  well  as  average  heat  transfer  coefficients  with  the  experimental  data  in 
Chapter  3.3.  This  value  is  therefore  used  for  all  of  subsequent  calculations. 


4.2.2.  Effect  of  Mass  Velocity 


Figure  4.5  shows  predicted  variations  of  several  film  parameters  with  quality  for 
different  FC-72  mass  velocities  at  a  constant  water  flow  rate  of  6  g/s.  Liquid  film  and 
vapor  core  Reynolds  numbers  are  defined,  respectively,  as 


'<2 

ii 

u  ' 

C4 

(4.50) 

and 

Re  =  P^. 

(4.51) 
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Figure  4.3  Effects  of  interfacial  dampening  term  on  (a)  eddy  momentum  diffusivity, 
(b)  velocity,  and  (c)  temperature  distributions  across  condensing  FC-72  film 
for  G  =  367  kg/  m2s,  and  mw  =  6  g/s. 
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(a) 


(b)  * 

Figure  4.4  Effects  of  interfacial  dampening  term  on  variations  of  (a)  liquid  film  thickness, 
and  (b)  local  heat  transfer  coefficient  with  quality  for  G  =  367  kg/  m2s  and  mK  =  6  g/s. 
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where  the  mean  liquid  film  velocity,  uf,  is  expressed  as 


'  Pr[W„,H,h-(Wah-2d\Hct-2S)\ 

Peterson  et  al.  (1997)  proposed  the  following  correlation  for  transition  from  laminar  to 
turbulent  film  flow  for  downward  flowing  films. 

Re fJr  =  4080[(l  +  0.3  lr*107  \\ .  1 8PrJ 87 - 0.23)]  ' ,  (4.53) 

where  the  nondimensional  interfacial  shear  stress  is  defined  as 


* 


PfM/gf 


(4.54) 


Based  on  this  correlation,  transition  to  turbulent  film  flow  is  estimated  to  occur  at 
very  low  film  Reynolds  numbers  when  the  values  of  interfacial  shear  and  liquid  Prandtl 
number  are  high.  Using  values  relevant  to  the  present  study,  with  r,.  =  15  Pa  (  r*  =  50) 
and  Tsal  =  60°C  (Pr/=  8.66),  Eq.  (4.54)  yields  a  transitional  Reynolds  number  of  Re/>  = 
25.  Therefore,  although  Re/  values  for  the  present  study  are  relatively  small  (less  than 
400),  strong  interfacial  shear  promotes  turbulence  film  flow  for  most  operating  conditions. 
This  fact  supports  assumption  (11)  used  in  the  development  of  the  present  model. 

Figures  4.5(a),  4.5(b)  and  4.5(c)  show  that  increasing  the  mass  velocity  of  FC-72 
increases  the  film  Reynolds  number,  vapor  core  Reynolds  number  and  interfacial  shear 
stress,  respectively.  The  slight  discontinuities  in  the  shear  stress  plots  are  caused  by  the 
different  ranges  of  the  effective  Reynolds  number  of  the  vapor  core,  Ret,  as  indicated  by 
Eqs.  4.21(a)  -  (c).  Figure  4.5(d)  shows  the  film  thickness  decreasing  with  increasing 
mass  velocity  because  of  the  increased  interfacial  shear.  Figure  4.5(e)  shows  how  the 
decreasing  film  thickness  with  increased  mass  velocity  increases  the  heat  transfer 
coefficient.  For  each  mass  velocity,  Fig.  4.5(e)  shows  that  the  heat  transfer  coefficient  is 
greatest  where  the  film  is  first  initiated  and  decreases  monotonically  downstream  as  the 
film  thickens  due  to  gradual  condensation  along  the  channel. 


4.2.3.  Validation  of  Model  Predictions 

Figure  4.6(a)  compares  pressure  drop  predictions  with  measured  pressure  drop 
data  for  different  water  mass  flow  rates  and  FC-72  mass  velocities  of  G  =  248,  306,  and 
367  kg/m2s.  This  figure  proves  that  the  model  accurately  captures  the  measured  pressure 
drop  in  both  magnitude  and  trend. 
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Figure  4.6(b)  compares  predictions  of  the  present  model  and  nine  prior  separated 
flow  correlations  against  the  present  pressure  drop  FC-72  data.  The  correlations  include 
six  specifically  developed  for  mini/micro-channels;  the  other  three  are  intended  for 
macro-channels.  All  three  macro-channel  correlations  significantly  over  predict  the  FC- 
72  micro-channel  data.  Among  the  mini/micro-channel  correlations,  those  by  Mishima 
and  Hibiki  (1996)  and  Nino  et  al.  (2005),  which  are  both  based  on  air-water  adiabatic 
flow,  show  fair  predictions.  Interestingly,  these  predictions  are  superior  to  those  of  the 
correlation  by  Tran  et  al.  (2000),  which  was  developed  specifically  for  flow  boiling  in 
small  channels.  This  supports  the  fact  that,  in  the  absence  of  droplet  entrainment  and 
deposition  effects,  micro-channel  condensing  flows  are  closer  in  flow  structure  (on  a 
local  basis)  to  adiabatic  than  boiling  micro-channel  flows  as  depicted  in  Fig.  3.8(a). 
Overall,  Fig.  4.6(b)  shows  the  excellent  predictive  capability  of  the  present  theoretical 
model,  evidenced  by  a  MAE  of  3.6%. 

Figure  4.7(a)  compares  predictions  of  the  present  model  and  average  FC-72  heat 
transfer  coefficient  data  for  different  FC-72  mass  velocities  and  water  flow  rates.  The 
model  does  capture  the  experimental  trends  relative  to  FC-72  mass  velocity  and  water 
mass  flow  rate.  Deviations  between  predicted  and  measure  values  may  be  related  to  film 
waviness,  which  is  more  prevalent  for  the  lower  mass  velocity.  Waviness  can  increase 
the  heat  transfer  coefficient  in  two  ways,  by  increasing  liquid-vapor  interfacial  area  and 
by  decreasing  mean  film  thickness  due  to  the  interfacial  waves  propagating  slightly  faster 
than  the  rest  of  the  liquid.  Another  reason  for  the  deviations  between  predicted  and 
measured  values  is  the  circumferential  non-uniformity  of  the  liquid  film,  with  more  liquid 
driven  towards  the  comers,  thinning  the  film  along  the  flat  walls  of  the  micro-channel. 
Despite  the  low  surface  tension  of  FC-72,  corner  effects  can  slightly  influence  the  mean 
film  thickness.  Another  source  of  predictive  error  may  be  minor  stratification  in  the 
horizontal  channels,  especially  at  low  mass  velocities. 

Figure  4.7(b)  compares  predictions  of  the  present  model  and  nine  prior  separated 
flow  correlations  against  the  present  average  heat  transfer  coefficient  FC-72  data.  The 
correlations  include  four  specifically  developed  for  mini/micro-channels;  the  other  five 
are  intended  for  macro-channels.  Interestingly,  excepting  the  correlation  by  Akers  et  al. 
(1958),  most  macro-channel  correlations  provide  good  predictions  of  the  present  data. 
With  a  MAE  of  4.7%,  the  recent  correlation  by  Kim  and  Mudawar  (in  review-a)  shows 
the  best  predictive  capability.  Nearly  as  accurate,  the  present  theoretical  model  predicts 
the  data  with  a  MAE  of  9.3%. 
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Figure  4.8  compares  predictions  of  the  present  model  and  prior  correlations  with 
the  present  local  heat  transfer  coefficient  data  for  three  FC-72  mass  velocities.  Notice 
that,  despite  their  low  MAE  in  predicting  the  average  heat  transfer  coefficient  data,  the 
correlations  of  Shah  (1979),  Cavallini  and  Zecchin  (1974),  and  Dobson  and  Chato  (1998), 
do  not  accurately  capture  the  variation  of  local  heat  transfer  coefficient  with  quality. 
Overall,  both  the  present  theoretical  model  and  the  recent  correlation  by  Kim  and 
Mudawar  (in  review-a)  show  good  accuracy  in  predicting  the  measured  trend  of  local 
heat  transfer  coefficient  with  quality. 
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Figure  4.6  (a)  Variations  of  predicted  and  measured  total  pressure  drops  with  water  mass 
flow  rate  for  different  FC-72  mass  velocities,  (b)  Comparison  of  measured  FC-72  total 
pressure  drop  data  with  predictions  of  present  model  and  previous  correlations. 
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Figure  4.7  (a)  Variations  of  predicted  and  measured  average  heat  transfer  coefficients 
with  water  mass  flow  rate  for  different  FC-72  mass  velocities,  (b)  Comparison  of 
measured  FC-72  average  condensation  heat  transfer  coefficient  data  with  predictions  of 
present  annular  model  and  previous  annular  condensation  heat  transfer  correlations. 
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Figure  4.8  Comparison  of  measured  FC-72  local  heat  transfer  coefficient  data  with 
predictions  of  present  model  and  previous  annular  condensation  heat  transfer  correlations 
for  mw  =  6  g/s  with  (a)  G  =  248  kg/  m2s,  (b)  G  =  306  kg/  m2s,  and  (c)  G  =  367  kg/  m2s. 
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CHAPTER  5.  CONCLUSIONS 


This  study  concerns  transport  phenomena  associated  with  condensation  in  micro¬ 
channels.  The  construction  of  the  condensation  module  and  flow  loop,  and  the 
experimental  methods  used  are  discussed.  High-speed  video  and  photomicrographic 
techniques  are  used  to  explore  and  help  categorize  two-phase  flow  regimes  associated 
with  condensation  of  FC-72  in  parallel  square  channels  with  Di,  =  1  mm.  The  two-phase 
pressure  drop  and  heat  transfer  characteristics  of  the  condensation  test  module  are 
examined.  An  assessment  of  the  accuracy  of  prior  models  and  correlations  in  predicting 
the  present  pressure  drop  as  well  as  the  heat  transfer  data  is  reported.  A  theoretical 
control-volume-based  model  tailored  to  turbulent  shear-driven  films  is  proposed,  and  the 
pressure  drop  and  heat  transfer  characteristics  of  annual  flow  condensation  in  rectangular 
micro-channels  with  three-sided  cooling  walls  are  examined. 


5.1.  Condensation  Pressure  Drop 

The  total  pressure  drop  measured  between  the  inlet  and  outlet  plenums  was 
compared  to  predictions  of  previous  models  based  on  homogenous  equilibrium  and 
separated  flow  models.  Key  findings  from  this  study  are  as  follows: 

(1)  Five  dominant  condensation  regimes,  smooth-annular,  wavy-annular, 
transition,  slug,  and  bubbly,  are  identified.  Transitions  between  these  regimes  are  in 
general  agreement  with  those  from  previous  studies.  Based  on  the  current  FC-72  flow 
visualization  data,  transition  lines  of  flow  regime  in  micro-channel  condensation  are 
proposed. 

(2)  The  total  pressure  drop  increases  with  increasing  mass  velocity  due  largely  to 
the  increased  interfacial  shear  stress  between  the  vapor  core  and  the  annular  liquid  film. 
Increasing  the  flow  rate  of  cooling  water  decreases  the  total  pressure  drop  slightly 
because  of  increased  flow  deceleration. 
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(3)  A  complete  model  was  constructed  for  total  pressure  drop  between  the  inlet 
and  outlet  plenums.  To  calculate  the  two-phase  frictional  pressure  drop  portion  of  the 
total  pressure  drop,  different  homogenous  and  separated  models  were  tested.  Although 
annular  flow  was  dominant  throughout  the  present  experiments,  most  homogeneous 
models  unexpectedly  provided  accurate  predictions  of  the  present  data.  All  separated 
flow  macro-channel  correlations  overpredicted  the  experimental  data,  some  with  very 
large  MAE.  Among  the  separated  flow  mini/micro-channels  correlations,  those  that  are 
based  on  adiabatic  and  condensing  two-phase  flows  provide  better  predictions  than  those 
based  on  flow  boiling.  This  can  be  attributed  to  fundamental  differences  in  local  annular 
flow  behavior  between  condensing  and  adiabatic  flows  on  one  hand  and  evaporating 
flows  on  the  other.  While  entrained  droplets  play  an  important  role  in  evaporating  flows, 
no  droplets  were  observed  in  the  present  condensation  experiments.  Overall,  weak 
predictions  of  separated  flow  correlations  highlight  the  need  of  further  study  and  new 
predictive  tools  specifically  tailored  to  micro-channel  condensing  flows. 


5.2.  Condensation  Heat  Transfer 

FC-72  condensation  heat  transfer  coefficient  data  were  compared  to  predictions  of 
previous  correlations,  and  a  new  correlation  was  developed  for  the  annular  regime.  The 
accuracy  of  this  correlation  was  ascertained  against  both  the  present  FC-72  data  and  a 
large  mini/micro-channel  database  from  eight  previous  sources.  Key  findings  from  this 
study  are  as  follows: 

(1)  The  FC-72  condensation  heat  transfer  coefficient  is  highest  at  the  point  of 
commencement  of  the  smooth-annular  regime,  where  the  annual  liquid  film  is  thinnest. 
The  heat  transfer  coefficient  decreases  along  the  micro-channel  because  of  gradual 
thickening  of  the  liquid  film  and  eventual  collapse  of  the  annular  regime.  For  high  mass 
velocities,  the  variation  of  the  heat  transfer  coefficient  is  flatter  for  intermediate  quality 
values  corresponding  to  the  wavy-annular  regime,  where  heat  transfer  enhancement  by 
the  interfacial  waves  partially  compensates  for  the  film  thinning. 

(2)  The  FC-72  condensation  heat  transfer  coefficient  decreases  with  increasing 
water  flow  rate  because  of  the  increased  rate  of  condensation  and  resulting  thickening  of 
the  annular  liquid  film.  Flowever,  the  condensation  heat  transfer  coefficient  is  less 
sensitive  to  variations  in  water  mass  flow  rate  than  in  FC-72  mass  velocity. 
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(3)  In  general,  prior  annular  condensation  heat  transfer  correlations  intended  for 
macro-channels  provide  better  predictions  of  the  present  average  FC-72  condensation 
heat  transfer  coefficient  than  more  recent  correlations  intended  specifically  for 
mini/micro-channels. 

(4)  A  new  condensation  heat  transfer  coefficient  correlation  is  proposed  for 
annular  condensation  heat  transfer  in  mini/micro-channels.  This  correlation  accurately 
predicts  the  present  FC-72  average  heat  transfer  data  with  a  MAE  of  9.0%. 

(5)  To  further  assess  the  predictive  accuracy  of  the  new  correlation,  a  database  of 
923  data  points  for  mini/micro-channels  (including  multi-channels  and  single  channels) 
was  amassed  from  eight  previous  sources,  of  which  713  data  points  correspond  to  annular 
flow.  The  new  correlation  shows  excellent  predictive  capability  based  on  the  713  annular 
data,  evidenced  by  an  overall  MAE  of  18.3%. 


5.3.  Theoretical  Model  for  Annular  Flow  Condensation 

A  theoretical  control-volume-based  model  is  developed,  which  is  based  on  the 
assumptions  of  smooth  interface  between  the  annular  liquid  film  and  vapor  core,  and 
uniform  film  thickness  around  the  channel’s  circumference.  The  model  predictions  are 
compared  with  experimental  micro-channel  pressure  drop  and  condensation  heat  transfer 
data  for  FC-72.  This  study  also  explored  the  accuracy  of  prior  separated  flow 
correlations  in  predicting  the  same  data.  Key  conclusions  from  this  study  are  as  follows. 

(1)  Unlike  flow  boiling,  droplet  entrainment  and  deposition  effects  are 
nonexistent  in  micro-channel  annular  condensing  flows. 

(2)  Annular  condensation  can  be  divided  into  two  separate  regions,  an  upstream 
smooth-annular  region  and  a  downstream  wavy-annular  region.  While  the  annular  film 
in  the  wavy-annular  regime  is  marred  by  interfacial  waves,  these  waves  are  fairly  small  in 
profile  and  do  not  evolve  into  large  waves,  nor  do  they  appear  to  have  an  appreciable 
influence  on  the  film’s  transport  behavior. 

(3)  For  shear-driven  films,  transition  from  laminar  to  turbulent  film  flow  may 
occur  at  unusually  small  film  Reynolds  numbers.  Turbulent  effects  must  therefore  be 
incorporated  when  modeling  the  transport  behavior  of  annular  condensing  films. 
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(4)  Interfacial  dampening  of  turbulent  eddies  is  accurately  accounted  for  in  the 
new  model  with  the  aid  of  a  new  eddy  diffusivity  model  specifically  tailored  to  turbulent 
shear-drive  films. 

(5)  While  some  of  the  previous  separated  flow  correlations  do  provide  accurate 
predictions  of  the  average  heat  transfer  coefficient,  they  fail  to  predict  axial  variations  of 
the  local  condensation  heat  transfer  coefficient. 

(6)  The  present  model  accurately  predicts  the  pressure  drop  and  heat  transfer 
coefficient  data  in  both  magnitude  and  trend,  evidenced  by  mean  absolute  error  values  of 
3.6%  and  9.3%,  respectively. 
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